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Abstract 

We present a detailed review of scaling behavior in the magnetically 
underdoped cuprate superconductors (hole dopings less than 0.20) and 
show that it reflects the presence of two coupled components throughout 
this doping regime: a non-Landau Fermi liquid and a spin liquid whose 
behavior maps onto the theoretical Monte Carlo calculations of the 2D 
Heisenberg model of localized Cu spins for most of its temperature do- 
main. We use this mapping to extract the doping dependence of the 
strength, f{x) of the spin liquid component and the effective interaction, 
Jeff{x) between the remnant localized spins that compose it; we find both 
decrease linearly with x as the doping level increases. We discuss the phys- 
ical origin of pseudogap behavior and conclude that it is consistent with 
scenarios in which the both the large energy gaps found in the normal 
state and their subsequent superconductivity are brought about by the 
coupling between the Fermi liquid quasiparticles and the spin liquid exci- 
tations, and that differences in this coupling between the 1-2-3 and 2-1-4 
materials can explain the measured differences in their superconducting 
transition temperatures and other properties. 
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1 Introduction 

Explaining the anomalous normal state properties of the so-called pseudogap 
regime of the underdoped cuprate superconductors is widely regarded as an es- 
sential step toward understanding the basic physics of these materials and un- 
locking the mechanism of their superconductivity[l]. Perhaps the most striking 
aspect of these is the universal, or scaling, behavior, first identified in mea- 
surements of their temperature-dependent uniform magnetic susceptibility [2], 
and since found in Knight shift, transport, and entropy measurements. In the 
present article we present a detailed review of scaling behavior in the underdoped 
cuprates that extends previous analyses of its manifestations in both static and 
low frequency dynamic behavior as well as that seen in inelastic neutron scat- 
tering (INS) experiments. Our review updates our earlier analysis [3] and the 
results presented in Norman et al.[l], and complements the recent review of gap 
behavior presented in Hiifner et al. [4] . 

We find that from zero hole doping until planar doping levels of ^ 0.2 are 
reached, the scaling behavior seen by probes of magnetic behavior reflects the 
presence of a spin liquid whose behavior maps onto the theoretical Monte Carlo 
calculations of the 2D Heisenberg model of localized Cu spins[5] for most of 
its temperature domain. We use this to extract the doping dependence of the 
strength, f{x) of the spin liquid component and the effective interaction. Jef f{x) 
between the remnant localized spins that compose it; for x < 0.18, f{x) = 
1 — [a;/0.20] for both the 2-1-4 and 1-2-3 materials, while, to first approximation, 
Jf,ff{x) = Jf{x), where J is their interaction at zero doping level. A careful 
analysis of the NMR experiments on both classes of materials makes it possible 
to identify a quantum critical point at a doping level, x = 0.05 that represents 
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a phase transition from short range to long range order in the spin hquid. It 
leads to dynamic z = I scaling behavior for a wide range of doping levels 
that extends up to T™, the temperature at which the static susceptibility is 
maximum, corresponding to an antiferromagnetic correlation length of order 
unity. We find that the extent of this scaling behavior is different for the 2-1-4 
and 1-2-3 materials: for the former it persists down to temperatures of order 
the superconducting transition temperature; for the latter it cuts-off at T*, a 
temperature that it considerably greater than Tc over most of the doping range. 

In addition to the spin liquid, whose properties dominate the low frequency 
magnetic response, bulk susceptibility measurements reveal the presence of a 
second component, a Fermi liquid that makes a temperature independent, but 
doping dependent contribution to this quantity for temperatures greater than 
T* and doping levels of 0.05 upwards. We present a simple interpretation of 
the two fluid description of these coupled liquids [3] in terms of the incomplete 
hybridization of the Cu d and O p bands; the spin liquid corresponds to the 
unhybridizcd d~ d component, while the Fermi liquid has a large Fermi surface 
as a result of the d — p hybridization. We derive the strength of the Fermi liquid 
component, which goes as [1 — f{x)] and so is proportional to x, and show how 
the presence of the spin liquid is incompatible with the single band Hubbard 
and Zhang- Rice approximations. 

We conclude that experiment has now provided the answer to the question of 
the physical mechanism responsible for the remarkable pseudogap behavior seen 
in the underdoped 1-2-3 materials {x < 0.20, say). When the present analysis is 
combined with the recent ARPES experiments [6] and the STM measurements 
of the Davis[7] and Yazdani[8] groups, a simple physical picture emerges. In 
the "normal state", for temperatures above T* ~ T^/3, one has two quasi- 
independent components: a spin liquid of localized Cu spins described by the 
2D Heisenbcrg model, whose strength and effective interaction become weaker as 
the doping level increases; and a (non-Landau) Fermi liquid with a large Fermi 
surface whose strength increases with doping and whose transport properties 
are determined primarily by its coupling to the spin liquid. At T* the system 
makes a transition to a remarkable new quantum state of matter: a state that 
possesses a single d-wave like gap, with a maximum gap value of order 4T* , 
that only becomes superconducting at the typically much lower superconducting 
transition temperature, Tc- The physical mechanism for the transition at T* 
(and subsequently at Tc) in the 1-2-3 materials is magnetic because the scale of 
T* and the gap is set by the effective interaction between the localized spins in 
the spin liquid. Matters are somewhat different for the 2-1-4 materials and we 
speculate as to why this is the case. 

Our review is organized as follows. In Section 2 we review the literature on 
experimental measurements and corresponding analyses that indicate univer- 
sal scaling behavior. In Section 3 we introduce the phenomenological two-fluid 
model and use it to analyze existing magnetic and thermodynamic measure- 
ments on the bulk spin susceptibility, the entropy, and the spin fluctuation 
spectrum revealed in nuclear magnetic resonance and inelastic neutron scatter- 
ing experiments. In Section 4 we present our conclusions concerning the inter- 
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action between the Fermi liquid quasiparticles and the spin liquid excitations, 
discuss the similarities and differences between the 2-1-4 and 1-2-3 materials, 
and consider the constraints imposed by experiment on microscopic theories of 
their high-tcnipcraturc superconductivity. We present our conclusions in Sec- 
tion 5. 

2 An overview of experiments suggesting uni- 
versal behavior 

The observation of scaling in the cuprates is not new [9]. Not long after the 
discovery of the cuprate superconductors [10], universal behavior was identified 
in the magnetic properties of the 2-1-4 materials by Johnston[2] through an 
analysis of his measurements of the bulk spin susceptibility; his analysis was 
later confirmed by Nakano et aZ.[ll] and Oda et al.[\2]. The Johnston-Nakano 
scaling analysis was subsequently extended by Wuyts et aZ.[13] to the Knight 
shift measurements of AUoul et aZ.[14] in the 1-2-3 family; more recently it 
has been shown to be applicable to the 1-2-4 and several other members of 
the 1-2-3 family by Curro et aL[15] and the authors[3]. A number of other 
experiments also indicate scaling and data collapse. These include electronic 
heat capacity measurements[16, 17] for which an analysis of the magnetic en- 
tropy found scaling behavior, quantum critical (QC) scaling behavior in NMR 
copper nuclear spin-lattice relaxation rates [18], limits on the T-linear behav- 
ior of resistivity[19, 13], scaling of Hall resistivity[20, 21, 13, 22], w/T scaling 
in inelastic neutron scattering[23, 24, 25, 26], and finite-size scaling in the in- 
sulating cuprates [27]. Recently, scaling behavior has been discovered in the 
doping-dcpcndcncc of ARPES (angle-resolved photoemission) and STM exper- 
iments on the 2-2-1-2 members of the 1-2-3 family; similar characteristic tem- 
peratures set the scale for the appearance of Fermi arcs [6], and "normal state" 
gap behavior [28]. 

In this section we review the above experiments and find that the character- 
istic scaling temperature first identified by Johnston, the temperature, T"^{x), 
at which the temperature and doping dependent bulk magnetic susceptibility 
reaches its maximum value, provides the common thread that links these to- 
gether. 

2.1 Direct measurements of the magnetic susceptibility 

The scaling behavior of the temperature dependent bulk spin susceptibility, 
x(T), was discovered empirically by Johnston[2], who showed that an excellent 
collapse (Fig.l) of his experimental data on five samples of h&2-z^T^zGnOi-y, 
in which the doping level x = z — 2y ranged from a; = to a; = 0.2, could be 
obtained if x(a;,T) had the following scaling form: 

x(x,T) = x^{x) + [x™(a;) ~ xo{x))]F {T / {x)) , (1) 
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where x is the doping level, Xo(a;) is a doping-dependent, temperature inde- 
pendent term, x™'{^) ~ x{x,T™-{x)) is the maximum value of x for a given 
doping level, and T"^{x) is the doping-dependent temperature at which x(r) is 
maximum. Johnston concluded that the scaling parameter T™ depends only on 
the hole doping level, x, in the plane and that the scaling function F{T/T"^) 
is the same as that calculated for the 2D Heisenberg model in its spin liquid 
regime (i.e. at temperatures above the Neel ordering temperature). In this 
model r™ = 0.93J, where J is the nearest neighbor exchange coupling be- 
tween localized Cu spins. The temperature independent Xo{x) was assumed 
to include x-independent core and Van Vleck contributions to Xi an x- 
dependent Fermi liquid contribution that grew somewhat slower than linearly 
with increased doping x. 

To explain the scaling behavior of Eq.(l), Johnston[2] introduced a doping- 
dependent 2D Heisenberg exchange constant, J{x), and the ratio 

where X^ici^) is the result obtained for the 2D Heisenberg model for a given 
J{x). He found that R{x) gradually decreases from R{x = 0) = 1 to R{x = 
0.2) ~ (Fig.2). 

Oda et al.[12] confirmed Johnston's scaling law for the bulk spin suscep- 
tibility in Sr- and Ba-dopcd La2Cu04, and extended its applicability to the 
2212 family of cuprates, Bi2Sr3_j,Caj,_2;Cu208. The experimental results for 
both families displayed excellent data collapse to the theoretical 2D Heisenberg 
curve. Like Johnston, Oda et al. found that the weight of the Heisenberg-like 
contribution decreases with hole doping. To account for this decrease, Oda et 
al. introduced an effective magnetic moment. 



Me// = gy/s{s + l). (3) 
g is defined in terms of T™{x) and x™{^) as 



They found that both T'^{x) and g{x) decrease linearly with doping, 

<?(a;) = 2.2(1 -4.5x) ocT"(x). (5) 

The magnitude of g{x) is that expected from the sum rule for a homogeneous 
material, 

r[X[^,q,x) - xo{x)) = g [x) , (6) 



27r y (27r?i) 

but we caution the reader that homogeneity may not be present in the under- 
doped regime. 
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Figure 1: Cu^+ sublattice susceptibility X^'^iT)/x^ax vs T/T"^''^ for five dif- 
ferent samples of La2-zSr2;Cu04_y with effective doping levels x ^ z — 2y = 0.0, 
X = 0.05, X ~ 0.1, X = 0.13, and x = 0.2, compared with theoretical 2D Heisen- 
berg calculations and the Weiss molecular-field (MFT) prediction. Reproduced 
with permission from [2]. 
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Figure 2: Effective intralayer Cu-Cu exchange constant J and the ratio R = 
Xmax I ')druix^'^ ^o\c doping p = X. Tc{p) is shown for comparison. The sohd 
curves are guides to the eye. Reproduced with permission from [2]. 



Both the doping dependence of the x-dependent Fermi Uquid part and the 
temperature dependent spin hquid component were found by Oda et al. to be 
in good agreement with ear her results of Johnston [2]; hke Johnston they found 
that the temperature-dependent scaUng part of the bulk spin susceptibility dis- 
appears at some critical doping value; according to Eq.(5), that is a; ~ 0.22. 

Another early study of bulk susceptibility in La2_a;Sra:Cu04 was performed 
by Yoshizaki et al. [29], who give a plot of T"^{x) ~ 1000i^(l - 4.5a;), without 
attempting Johnston scaling. Since lOOOA' does not match the value of the 
exchange coupling J in the insulator, Yoshizaki et al. found a jump in T™{x) 
at the MI boundary. 

Further confirmation of Johnston scaling of the form Eq.(l) for the bulk 
spin susceptibility in the 2-1-4 family was obtained by Nakano et al. [11] (Fig. 
3), who demonstrated an excellent data collapse for a number of samples of 
La2-xSra;Cu04, both in the underdoped region and that close to and beyond 
optimal doping. Nakano et al. arrived at their scaling law by assuming the pres- 
ence of additional temperature-dependent terms: an impurity Curie term C/T 
and a linear term B{T — To) in the underdoped and overdoped regimes. They 
did not attempt to fit the 2D Heisenberg model calculations; an empirical scal- 
ing data collapse was constructed instead, with results that were in agreement 
with Oda et a/. [12]. 

An alternative form of scaling for the bulk spin susceptibility was proposed 
by Levin and Quadcr[30]. Similar to the studies reviewed above, they suggested 
separation of the bulk spin susceptibility in two components, the temperature- 
independent Fermi liquid component and the scaling component. However, in 
their model the scaling component originates from the contribution of a separate 
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Figure 3: X^{T)/Xmax versus T/T^ax for the superconducting samples of 
La2-a;Sra;Cu04. The arrow shows the value of x^{T = ^)/Xmax predicted by 
the 2D Heisenbcrg model calculations. The solid line is the universal curve. The 
inset shows that the universal curve decreases more rapidly for T <^ Tmax than 
that obtained by Johnston[2] (dashed line). Reproduced with permission from 
[11] 
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itinerant band. Levin and Quader[30] took the 2D density of states in the 
following form: 

iy{e) ^ ly^ + i^r^ei-e), (7) 

where the index ^ corresponds to a large hole band, which produces the usual 
Fermi liquid term in the bulk spin susceptibility. The chemical potential /x 
lies very close to the top of the second band 77 and enters that band at a; = 
xo- The contribution of the band rj to the bulk spin susceptibility becomes 
temperature-dependent due to thermally activated carriers and is a universal 
function Xni^ix — Xo)/T), where F cx is a Fermi liquid parameter. The pa- 
rameter T{x — xq) is thus an analog of T™-{x). The two-band model produced a 
reasonably successful data fit for the bulk susceptibility in TlSr2(Lui_j;Ca2;)Cu20y 
for both the underdoped and the overdoped regimes. 

To summarize, the various independent measurements of the bulk suscepti- 
bility in the 2-1-4 and 1-2-3 materials can now be seen to be consistent with 
one another and with the picture first set forth by Johnston: that in the under- 
doped regime, for dopings between x ~ 0.06 and ^ 0.22, and T > r™/3, one 
has two independent contributions to the bulk spin susceptibility. One comes 
from a 2D Heisenberg spin liquid with a doping-dependent effective interac- 
tion, Jeff{x) ~ T"^{x); the second represents a Fermi liquid contribution whose 
strength increases as the doping level is increased. As we shall see, the fall- 
off of the rescaled spin susceptibility below the Heisenberg spin liquid value at 
T ~ T™/3 appears to be a universal property of the spin liquid in the cuprates. 

2.2 Measurements of the bulk susceptibility using the NMR 
Knight shift 

The Knight shift seen in NMR experiments is a measurement of the bulk spin 
susceptibility at a particular nuclear site[31], so that measurements of the Knight 
shift for different nuclei serve to supplement direct measurements of the bulk sus- 
ceptibility. Most of the early analysis of the NMR data used a single-component 
Mila-Rice-Shastry (MRS) description [32] for which the justification was the ob- 
servation by AUoul et aL[14] and Takigawa et aL[33] that the ^^Cu, ^^O and 
^^Y Knight shifts in YBa2Cu307 and YBa2Cu3 06.63 have the same anomalous 
temperature dependence (Fig. 4). MRS proposed a hyperfine Hamiltonian that 
described the coupling of a single magnetic component formed by the system 
of planar Cu^"*" spins and holes mainly residing on the planar copper sites to 
the various nuclei. Most earlier Knight shift experiments[33, 34] confirmed this 
one-component Zhang- Rice [35] singlet picture, which is basically correct for the 
parent insulator. However, as first noted by Walstedt et al.[36] in connection 
with spin-lattice relaxation rate measurements, and discussed in detail below, 
the single component description turns out to lead to a number of contradictions 
in the doped materials, and requires modification. 

Just as was the case for the bulk spin susceptibility, the temperature-dependent 
Knight shift data for different doping levels and different families of cuprates 
displays Heisenberg model-type scaling. The ^^Y Knight shift data from Alloul 
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Figure 4: The experimental planar Cu and O Knight shifts in YBa2Cu306.63 
plotted vs the temperature T . The values for the static planar susceptibility 
Xo(T')/Mb are given on the right-hand scale. Different Knight shifts and the 
static spin susceptibility have the same temperature dependence. Reproduced 
with permission from [34]. 



et al. [14] for different doping levels in the 1-2-3 family has been scaled to a single 
curve with very good data collapse by Wuyts et aL[13], although no comparison 
to the 2D Hcisenberg model or the Nakano et a/. [11] bulk spin susceptibility 
scaling curve was provided. Recently Curro et a/. [15] and the authors[3] have 
shown that the two are consistent; as may be seen in Fig. 5, the NMR Knight 
shift follows very well the Johnston-Nakano bulk susceptibility scaling form and 
the Heisenberg model for the 1-2-4 member of the 1-2-3 family. Spin liquid 
scaling behavior determined by T"^{x) has thus been shown to be universal in 
the cuprates. 

What then is the direct experimental evidence from Knight shift measure- 
ments for the presence of two distinct components? One material for which a 
contradiction with the simple one component model was found is YBa2Cu408, 
where a more precise measurement [38] of the copper Knight shift data in mag- 
netic fields parallel to the c-axis, ^^K\\, found that it displays a rather un- 
usual temperature dependence, one that is different from that seen for the bulk 
spin liquid susceptibility. A second example comes from quite recent measure- 
ments by Haase et aZ.[39] of the planar and apical oxygen Knight shifts for 
Lai.85Sro.i5Cu04. As we discuss in a later section, these indicate that while 
the temperature-dependent part of the Knight shift on different nuclei is indeed 
the same, there is a temperature-independent part of the Knight shift above the 
superconducting temperature, T^, that is different for the planar copper and the 
planar and apical oxygen nuclei. The key signature of this effect, the deviation 
of various Knight shifts from the same temperature dependence at tempera- 
tures below Tc, is clearly visible in Fig. 4. These contradictions suggest that the 
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Figure 5: The spin shift in YBa2Cu408 versus T/T^"-, showing the apphcabihty 
of Johnston/Nakano scahng and the relevance of the 2D Heisenberg modcl[15]. 
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ceptibility of the spin-1/2 2D Heisenberg antifcrromagnct [5] and the data points 
are taken from Ref.[37]. The arrow points to the zero-temperature prediction 
for the 2D Heisenberg antiferromagnet. Reproduced with permission from [15]. 
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original MRS hyperfine Hamiltonian has to be modified to include these new 
effects, and that the one-component spin dynamics and the Zhang-Rice singlet 
picture [35] fails at moderate doping levels [36], a topic to which we return in 
Section 3. 



2.3 Spin-lattice relaxation rates 

In the MRS Hamiltonian that is described in detail in the following section, wave 
vector dependent form factors arise from the presence of a transferred hyperfine 
interaction between the Cu spins and the probe nucleus. For probe nuclei other 
than copper, these vanish at the commensurate wave vector Q = (tTjTt), so 
that, for example, an in-plane oxygen nucleus will feel little of a spin response 
that is peaked at the commensurate wave vector. (That this should be the 
case is obvious if one recalls that such oxygen are located midway between 
copper nuclei, and for an antiferromagnetic array of copper spins, the nearest 
neighbor spins would cancel one another out in their influence on the oxygen 
site.) The striking difference of the temperature dependence and magnitude of 
the spin-lattice relaxation on the ^^Cu, ^^O, and ^^Y nuclei led Millis, Monien, 
and Pines (MMP)[40] to the conclusion that a localized or nearly localized spin 
component of the dynamic magnetic response function must be strongly peaked 
at the commensurate wave vector, as might be expected if one were close to an 
antiferromagnetic instability. 

Indeed, closer study shows that even a slight deviation from commcnsurabil- 
ity within the MMP approach based on the MRS one-component model will have 
a significant impact on oxygen relaxation rates that is not seen cxperimentally[40] . 
Thus, an incommensurate peak structure for x"(q, w), such as has been inferred 
from inelastic neutron scattering (INS) experiments[41, 42] on the 2-1-4 mate- 
rials, is inconsistent with this approach. One way to get around this difficulty 
is to introduce additional transferred hyperfine interactions [43]; a second way, 
proposed by Shchter[44, 45] , is to note that unlike NMR, INS is a global probe 
of spin excitations, so that a suitable domain structure (regions of commensu- 
rate near-antiferromagnetic behavior, separated by domain walls) would give 
rise to the apparent incommensuration inferred from the INS experiments. We 
adopt this explanation in what follows. 

The spin fluctuation response function proposed by MMP was that appro- 
priate to any spin liquid near a commensurate antiferromagnetic instability: 

XSL{q,u;)^ ^ , (8) 

where the peak susceptibihty takes the form, 

XQ = af, (9) 

with ^ as the magnetic correlation length, and a{x) as a temperature- independent 
constant. The copper NMR ^^Ti and ^^T2g relaxation rates provide a direct 
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measure of the strength and character of the spin hquid response function, mea- 
suring as they do the momentum-integrated imaginary and real part of the spin 
response function, xsi(q, w)[46]. One finds: 

' (10) 



63Ti iUSF 
1 



63 



T2G 



OC (11) 



Imai et al.[l8\ in an early scaling analyses (Fig. 6) of the copper relax- 
ation rates for the 2-1-4 family demonstrated experimentally the existence of 
a universal high temperature limit for ^'^Ti that is temperature- and doping- 
independcnt: 

^^T^) = 

which according to Eq.(lO) means that at high temperatures the spin fluctuation 
energy, losf must be proportional to T. This universal high-temperature behav- 
ior of ^^Ti was explained in terms of the QC (Quantum Critical) [46, 9] behavior 
seen in the sigma-model[47, 48] that, strictly speaking, is only applicable to the 
parent insulating compound. Indeed the empirical finding of a temperature- 
and doping independent ^^Ti in the 2-1-4 family of materials extends to tem- 
peratures much higher than those at which such a theoretical explanation would 
apply. 

Two other high-temperature forms for the relaxation rates based on 2 = 1 
and z — 2 dynamical scaling were therefore suggested [46, 9] on the basis of 
the non-linear sigma modcl[47, 48] and general scaling arguments, and verified 
experimentally. In z = 1 scaling, ujsf and ^ are related by 

LOSF OC A = - , (13) 

so the presence of z = 1 QC scaling in the underdoped materials leads to the 
simple results: 

63^ OC USF^ OC C (14) 

On the other hand, for z = 2 (mean field) scaling one has 



r , , 

WSF = , (15) 



so that if it is present one finds 

63^^ 2^ 



637^2 



cx wsFr = r. (16) 



The NMR experimental results of Curro et aL[37] show that both forms of 
scaling are present in the material, YBa2Cu408. At high temperatures, for 
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T > T™ ~ 500K, Eq.(16) is valid and one has mean field behavior, while below 
T"' the spin spectrum displays z = 1 QC behavior down to a temperature 
~ 150K, where the Ti measurements suggest that a gap opens up in the spin 
liquid spectrum. We return to this finding below. 

The relaxation rates for other nuclei, such as oxygen or yttrium, also display 
anomalous (i.e. non-Korringa) behavior, although in much milder form. A 
modified Korringa-type scaling for oxygen ^^Ti was suggested by MMP: 

"rirxo(T) = const. (17) 

Additional scaling forms have been proposed for the copper relaxation rates 
at lower temperatures. In particular, the authors[3] have observed recently 
that ^^TiT reaches its universal high-tcmperature behavior with a different x- 
dependent offset (Fig. 7), 

^^TiT = BT + A{x), (18) 

where the constant B is universal, while A{x) changes linearly with doping, and 
its variation suggests the existence of a QC point in the spin liquid at xq — 0.05, 
where A{xo) — 0. Eq.(18) shows excellent data collapse for the 2-1-4 and the 
1-2-3 families. 

An alternative form of scaling has been suggested to apply at moderate to 
low temperatures, T < 300A', by Gor'kov and Teitel'baum (GT)[50] who found 
excellent data collapse (Fig. 8) for the following scaling form: 

(19) 



63Ti ^'^Ti{x) 63f^(r) ' 

where Ti(T) is a universal function for all high-Tc materials, while 637^"|(j;) 
varies with x. The suggested scaling form follows from the decomposition of re- 
laxation rate into two different processes. The x-dependcncc of l/^'^Ti(a;)(Fig.9) 
looks rather unusual[50], since the proposed empirical disorder-driven relaxation 
rate, l/^^Ti{x) decreases with increased doping x. The proposed GT scaling 
is assumed to be the result of intrinsic phase separation, and the pseudogap 
temperature T*{x) is proposed as a measure of the onset of such behavior. 

2.4 Finite size effects 

If scaling is present, it is natural to look for the possible presence of finite size 
effects, as Cho et al.[21] have done in the lightly doped La2-xSr2,Cu04 material. 
They found very good agreement with that expected for phase separation and 
domain formation on the insulating side of the Mott transition. However, there 
was no direct observation of the domain sizes of the different phases. According 
to finite-size scaling theory [51], the Neel temperature for a domain of size L 
takes the form, 

1 - « (20) 

1n{L = oo) 
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Figure 8: Temperature dependence of for YBa2Cu3 06.63 overlaid with 

that for YBa2Cu408 and Lai.8.5Sro.i5Cu04. In the lower right corner: the offset 
values for different compounds. Reproduced with permission from [50] 




Figure 9: The offset if^^Tiix) vs Sr content x for LSCO (relative to that for 
YBa2Cu408). The line is a visual guide. Reproduced with permission from [50]. 
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where v = 1/2 m. mean field theory. The experimental data can also be fit well 
by the expression 

with n = 1.90 ± 0.20 ~ 2. Cho et al. conclude that 

L{x) oc - , (22) 

X 

which means that the width of the domain wall is x independent. They found 
that in the antiferromagnetic regime the scaling law takes the form: 

x{x,T)=x{F{x){T ~Tn{x))), F{x)cx- . (23) 



' Co(T«(3n)) 

Here S,o{T) is the correlation length in the pure Hcisenbcrg model, 

^o(T) = 0.276a exp(27rp,/fcsr) 0.276a exp {J /ksT). (24) 



2.5 Inelastic neutron scattering 

Inelastic neutron scattering (INS) experiments enable one to explore the extent 
to which 2D Heisenberg model captures the momentum dependence and be- 
havior at higher frequencies of the spin liquid as its properties are altered by 
doping. There is by now a vast body of literature that includes several recent 
reviews [52, 53, 54]. Our focus here is on the extent that the doped spin liquid 
continues to exhibit the quantum critical behavior inferred from the ultra low 
frequency dynamic properties measured in the NMR experiments described in 
the previous sections, and the ways in which departures from quantum critical 
behavior emerge as the frequency is increased into the multi-time range or one 
goes into the gapped normal or superconducting state. As we shall see, one of 
the most striking features that emerges with doping is the observation of peaks 
whose positions at lower frequencies reflect a doping-dependent incommensura- 
tion or discommensuration that indicates dynamic stripe formation; a second is 
the appearance of resonances and spin gaps in the normal state. Throughout 
this section we will be concerned, as we were in our discussion of NMR experi- 
ments, with possible universal behavior that is common to the 1-2-3 and 2-1-4 
families. 



2.5.1 uj/T scaling for the local spin susceptibility 

Perhaps the most direct evidence of quantum critical scaling was provided by 
the INS measurements of the local (integrated) dynamic spin susceptibility. 
The remarkable results obtained in insulating and metallic low doped samples 
indicate the presence of the to/T scaling expected in the vicinity of a z = 1 
quantum critical point. 
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Figure 10: Normalized integrated spin susceptibility as a function of the scaling 
variable lo/T. The solid line is the function 2/7r tan~^ [oiw/T + a3(ijj/r)'^] with 
ai = 0.43 and 03 = 10.5. Reproduced with permission from [23]. 



Among early neutron scattering experiments that exhibit some form of scal- 
ing behavior those of particular interest are on lightly doped YBa2Cu3 06+a: by 
Birgeneau et al. [24], lightly doped La2-a;Sr^Cu04 by Keimer et al. [23], and 
YBa2Cu306.6 by Sternliebe< al. [25]. Keimer et aL[23] and Birgeneau et aZ.[24] 
fit their experimental results for the local (integrated) spin susceptibility to the 
expression 

X'/H^ j ^x"{q,u;)=C{u;)tan-\a^{co/T) + a,{cj/Tf + ---) (25) 

The data collapse is good (Fig. 10), but this form does not quite exhibit u/T 
scaling because of the T-independent amplitude C{uj). 

Sternlieb et al. [25] subsequently found that true uj/T scaling exists for 
the local spin susceptibility in the underdoped material YBa2Cu306.6- Their 
experimental results (Fig. 11) can be fit to the following simple scaling form: 

xU^,T) = Aix)^, (26) 

with a deviation from scaling that occurs at progressively lower temperatures, 
as the frequency to increases. 

In more recent experiments, Bao et al.[55\ studied experimentally the de- 
struction of 2D antiferromagnetic order in Li-doped La2Cu04. Since holes are 
loosely bound by Li impurities, the doped material remains an insulator even 
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Figure 1 1 : (top) The temperature dependence of the local dynamic spin suscep- 
tibility x"{'-^) for a fixed energy transfers converges to a universal function of 
T/lo at high temperatures. The arrows show the temperature above which the 
scaling behavior occurs, (bottom) The corresponding scattering function S{iu) 
vs T/io. Reproduced with permission from [25]. 
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Figure 12: (a) uj/T scaling is valid for La2Cuo.94Lio.06O4 in the high temperature 
QC regime The solid line is the scaling function, Eq.(28). (b) lo/T scaling 
becomes invalid in the low temperature regime, (c) A new scaling for the low 
temperature regime, with a constant energy scale Fo = ImeV. The solid line is 
the scaling function Eq.(30). Reproduced with permission from [55]. 



when long-range order is destroyed at x > 0.03. However, unlike the supercon- 
ducting La2_a;(Sr,Ba)a;Cu04, there are no additional complications due to the 
presence of mobile doped holes, making possible a cleaner measurement of the 
spin dynamics. Bao et al.[55] found a commensurate energy spectrum of spin 
excitations at a; > 0.03, where the long-range order is not present, and a charac- 
teristic quantum critical uj/T scaling for the local spin susceptibility (Fig. 12). 
The scaling function, however, is different from that obtained by Sternliebet 
a^.[25] In particular, Bao et aL[55] found significant deviations of the scaling 
function from linearity, with 

xl{u;,T)^X.f{^m, (27) 

and 

Below T = 50K this scaling changes to one with a constant energy scale Fo = 
ImeV, 

xULo,T) = x.9ito/To), (29) 

where 

.(.) = (30) 

Stock et al. [56] recently confirmed the low- frequency w /T scaling for x'l ('^) ^) 
in oxygen ordered ortho-II YBa2Cu306.5 superconductor found earlier in this 
material in the oxygen-disordered state by Birgeneau et al. [24] . Stock et al. fit 
their expression to the form Eq.(25), with only a linear term in uj/T present. 
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Figure 13: Temperature dependence of the inverse correlation length k{to,T) at 
various fixed energy transfers plotted against y^^aJ/fcep+T^. The solid line 
corresponds to a z = 1 QCP. The graph in the upper right shows how the peak 
response depends on fc = k{uj ~ 0,T). The inset in the upper left shows the 3D 
space defined by {uj,T) phase space probed by the x = 0.14 sample, and the 
solid circle represents the nearby QCP. Reproduced with permission from [26]. 



They find that the uj/T scaling breaks down at higher frequencies, lo > 20meV, 
since the amplitude C'(lu) in Eq.(25) becomes temperature-dependent. 

2.5.2 Temperature and frequency dependence of the correlation length 

Keimer et a/. [23] find evidence for finite size scaling, since a good fit to their 
experimental results can be obtained with the following expression for the cor- 
relation length: 

^-\x,T) = r\^,T^ 0) + = 0,T), (31) 

with ^{x,T 0) oc l/x (X where L is the size of the domain for finite size 
scaling. 

Eq.(31) was later generalized to finite frequencies u! by Aeppli et al.[26], who 
fit their results for the correlation length in Lai.g6Sro.i4Cu04 (Fig. 13) to the 
following form: 

r'{T, uj) = v/e-'(0)+c.2 + r2. (32) 
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2.5.3 Direct measurements of x"(g,ti;) 

Inelastic neutron scattering experiments provide a direct measurement of the 
spectrum of spin excitations, and thus of the spin wave velocity and exchange 
couplings in the antiferromagnetic insulator[57, 58, 59]. We shall see that the 
detailed measurements of the spectrum of the spin excitations in underdoped 
high-temperature superconductors that have recently been carried out for both 
the 2-l-4[60, 61, 62, 63] and the l-2-3[64, 65, 66] families provide evidence for 
its universality. 

Incommensuration/Discommensuration Early neutron measurements on 
both the 2-1-4 and the 1-2-3 families of high-tcmpcraturc superconductors re- 
vealed an incommensurate spectrum of spin excitations at low frequencies. 
Cheong et a/. [41] and Mason et a/. [42] discovered in their inelastic neutron scat- 
tering measurements on the 2-1-4 family of materials that the position of the 
low-energy spin fluctuation peak shifts from (tt, tt) in the parent insulating com- 
pound to (7r±27r(5, tt) and (tt, tt±2t:5) in the underdoped superconductor. They 
found that the incommensurability S in the underdoped regime grows approxi- 
mately linearly with increased hole concentration x, 6 x, at least for x < 0.12. 
For X > 0.12, the incommensurability S{x) saturates at a finite value. Incom- 
mensurate peaks in the spin response function x"{q,^) at low frequencies were 
also discovered in the bilayer system YBa2Cu306+i^[67]. As Mook et al.[67] 
have first shown, the incommensurate peaks in the bilayer system are located at 
the same positions as in the single-layer La2_xSra;Cu04 system at similar level 
of hole doping. Thus, they concluded that the band structure, which is very 
different in these two families of materials, plays a minor role in the spectrum of 
spin excitations at low frequencies. In what follows we will refer to these peaks 
as reflecting incommensurate behavior, although, as noted earlier, for these re- 
sults to be consistent with NMR, these should rather be regarded as reflecting 
discommensuration, or dynamic stripe order. 

Wakimoto et aL[68] and Matsuda et a/. [69] extended the measurements of 
x"{q,ijj) in La2_i;Sr3;Cu04 to the spin glass regime, 0.02 < x < 0.055. They 
found that while the low-energy spin excitations remain incommensurate with 
the same J ~ x as in the metallic phase, the positions of the incommensurate 
peaks arc rotated by 7r/4 relative to those in the metal, as shown in Fig. 14. 
Unlike in metallic phase, the diagonal stripe structure for the insulating spin 
glass phase is seen as Bragg peaks in elastic scattering measurements and thus 
is static at low temperatures [69]. Moreover, the incommensurability of the spin 
fluctuations disappears at higher frequencies or temperatures, above a certain 
energy threshold for the incommensurate structure, which we shall see is directly 
related to the spin gap. 

An interesting linear scaling relationship between the incommensurability 
d{x) and the superconducting temperature, T^^x), was suggested by Yamada et 
al.[70]. which, they find, holds extremely well for the 2-1-4 family of materials 
in the underdoped regime, Tc{x) — hv *2i4 S{x), with hv*2ii — 20meVA. Bal- 
atsky and Bourges[71] and Dai et aL[72] found that Yamada scaling is also 
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Figure 14: Hole concentration (x) dependence of the splitting of the incom- 
mensurate peaks (S) in La2-a;Sra;Cu04 in tetragonal reciprocal lattice units. 
Open circles indicate the data for the inelastic incommensurate peaks reported 
by Yamada et aL[70]. Filled circles and square are the data for elastic incom- 
mensurate peaks reported by Wakimoto et al.[68]. The broken and solid lines 
correspond to S = x. The insets show the configuration of incommensurate 
peaks in the insulating phase (diagonal stripe) and the superconducting phase 
(collinear stripe) Reproduced with permission from [69] . 



applicable to the 1-2-3 family of materials, with Tc{x) ~ fiv *i23 5{x) and 
t'*i23 ~ 36.6meT^A[72]. Since the maximum of Tc is different in different fami- 
lies of cuprates, while the incommensurability is universal and depends only on 
hole doping, is different for different families of materials. 

Stock et aZ.[56] considered the applicability of Yamada scaling to both 2-1-4 
and 1-2-3 families of high-temperature superconductors. They plotted the in- 
commensurability (5 as a function of Tc(a;)/r™°'^, a quantity proportional to the 
number of holes in the plane in the underdoped regime (see Fig. 15). Their plot 
thus confirms for the 2-1-4 materials the conclusion of Mook et aZ. [67] that the 
incommensurability of low-energy spin fluctuations in cuprate superconductors 
depends only on hole concentration in the plane, and not a particular material 
or band structure. 
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Figure 15: Reduced superconducting temperature Tc/Tcijnax) as it relates to 
the incommensurability S in YBa2Cu306+a: and La2-a:Sra:Cu04. Reproduced 
with permission from [56]. 
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High-energy spin excitations and the resonance peak Early INS studies 
of high-frequency spin fluctuations were limited due to poor resolution. Many of 
these measurements were initially done on the 1-2-3 family of materials, where 
the resonance peak[73, 74] enhances spin response at higher energies. The reso- 
nance peak in the 1-2-3 family of materials corresponds to a sharp enhancement 
of the intensity of commensurate spin excitations at high frequencies near the 
frequency, Eres{x), that depends on hole doping x. The observation then raised 
the question of whether the peak was a consequence of spin gap in the fermionic 
excitations or was an intrinsic property of the spin liquid in the vicinity of (tt, tt). 

At first the resonance was only observed in the superconducting state. Since 
the resonance peak was not observed in the inelastic neutron scattering exper- 
iments on the 2-1-4 family, this led to an explanation that this feature most 
likely arises as a result of the coupling of spin excitations to fermions with a 
d-wave gap or d-wave like gap in the normal state. Indeed, Dai et aZ.[74] found 
that the resonance peak first appears above Tc at T*{x), which decreases with 
increased hole doping, the resonance energy Eres{x) increases with increased 
doping and tracks the doping dependence of Tc{x), as shown in Fig. 16. 

However, more detailed neutron scattering studies revealed that the reso- 
nance corresponds to a special frequency at which high-energy spin excitation 
spectrum becomes commensurate. As we have seen above, the spin excitation 
spectrum depends only on hole doping, and not on a particular material, a 
conclusion that casts doubt on the quasiparticle gap explanation of the reso- 
nance. Inelastic neutron scattering studies near the resonance frequency were 
first done by Bourges et aL[75, 76] in the YBa2Cu306.5 material. Bourges et al. 
found that the commensurate resonance peak was broadening in momentum, 
both above and below the resonance frequency Eres- Based on their findings, 
Bourges ei al. suggested that the spin excitations disperse at high frequencies 
in a way that is similar to that expected for spin waves in the parent insulating 
compound. Aral et aL[77] later studied the momentum dispersion of spin ex- 
citations near the resonance peak in YBa2Cu306.7 in more detail. They found 
evidence for two modes that meet near the resonance energy, one that opens 
downwards and gives rise to incommensurate spin excitations at low energies, 
the other, a new mode that opens upwards in energy, giving rise to spin wave- 
like dispersive excitations at high energies. Arai et al. found that the two modes 
meet at the frequency to ~ AlmeV , slightly above the characteristic resonance 
frequency Eres = 36meV for their material. More recent INS experiments on 
YBa2Cu306.5[65], YBa2Cu306.6[64] and YBa2Cu306. 95(66], however, suggest 
that the high energy mode and the low energy mode meet at the resonance 
frequency. 

A typical energy spectrum of the spin excitations as measured by Stock et 
aZ. [65] in a detwinned ortho-II sample of YBa2Cu306.5 is shown in Fig. 17. As 
Mook et al.[78] had found earlier for a detwinned sample of YBa2Cu306.6j Stock 
et al. found two, rather than four incommensurate peaks below the resonance 
frequency Eres ~ 33meV^. This result strongly suggested a one-dimensional 
character for the spin fluctuations, consistent with the formation of dynamic 
stripes [79, 61]. The incommensurability 5 observed at low frequencies decreased 
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Figure 16: T* and the resonance energy Eresonance as a function of Tc in 
YBa2Cu306+2:. The open circles and open squares are temperatures at which 
dp(T)/dT reaches broad maximum. The open diamonds show the pseudogap 
temperature T* determined from NMR measurements [33]. The filled circles cor- 
respond to Tc and T*, where the resonance first appears in INS measurements. 
Filled squares are from Fong et al. [73] . Horizontal error bars are superconduct- 
ing transition widths. The solid lines are guides to the eye. Reproduced with 
permission from [74]. 
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Figure 17: The dispersion of acoustic and optic modes in YBa2Cu306.5 with 
respect to commensurate (tt, tt) position at T = 6K. The fiUed circles are from 
two Gaussian fits along the [100] and [010] directions above the resonance energy 
and from [100] direction below the resonance energy. The open circles represent 
the positions of the incommensurate peaks found in experiments conducted at 
Chalk River. The open squares represent the peak position that results from a 
polar average around (tt, tt) position. The solid lines schematically represent the 
dispersion of the insulating compound as measured by Hayden et al.[58\. The 
dashed lines are fits of the high-energy dispersion to linear spin wave theory for 
Lo > AOmeV. The horizontal bars show the q-width observed for three optic 
scans at constant energy. Reproduced with permission from [65] . 
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with increased energy transfer and disappeared at w = Eres- The spin excita- 
tions started to disperse outward again at high frequencies, lu > Eres- Stock 
et al.[65] found the dispersion for acoustic spin excitations to be isotropic and 
similar to the spin waves in the insulator, as shown in Fig. 18; they fit their 
results to a gapped spectrum, 

eaciq) - VAirTM^, (33) 

with Aqc — Eres — iSmeV and a spin wave velocity c ~ 365meVA^ which is 
close to the value c ~ AOOmeVA obtained from the slope of high energy spin 
excitations. The spin wave velocity is dramatically reduced from its value in 
the insulator, c ~ 650meyA[58]. 

Stock et al. also measured the width of the ring of high energy spin wave 
excitations and found that it increased linearly with frequency at low temper- 
atures. The low-energy spectral weight was dominated by the resonance peak 
at Eres = 33meV^, which had an asymmetric shape, with a quick drop-off of 
intensity for lo > Eres and a much slower reduction at oj < Eres- Stock et 
al.[56], however, estimated that the resonance represented only 3 % of the total 
spectral weight of spin excitations; most of its spectral weight would appear only 
at higher energies. Similar to spin waves in the insulating compound, x"l(w) 
approached a constant at high energies to > Eres, > Eres) — 5 — 

somewhat reduced from a similar value for the insulator, 7 — 8^'^/eV[65\. 

A very similar dispersion for spin excitations was found by Reznikei aZ.[66] 
in a twinned sample of YBa2Cu3 06.95- They observed four incommensurate 
peaks dispersing inward with increased frequency at frequencies below the reso- 
nance frequency Eres = 'ilmeV and becoming commensurate at the resonance 
frequency. The dispersion of spin excitations was cut off below the resonance 
frequency at a characteristic spin gap frequency Egap — SSmeV. Reznik et 
al. concluded that their data was consistent with the isotropic gapped spin 
wave-like excitations dispersing outward at frequencies above the resonance fre- 
quency. As had Stock et al., Reznik et al. found the q-integrated intensity to 
be approximately constant at high frequencies co > Eres, as expected for spin 
wave-like excitations. Because of resolution problems, Reznik et al. were not 
able to calibrate their data in absolute units, or estimate the spin wave velocity. 

Not all INS studies have found spin wave-like excitations with a reduced 
effective exchange coupling J at high energies. For example, Hay den et aZ.[64] 
performed measurements in a twinned sample of YBa2Cu306.6 with effective 
hole doping x ~ 0.1 close to the x = 1/8 stripe ordering instability. As had 
other groups, Haydcn et al. observed incommensurate spin excitations that 
dispersed inward up to resonance frequency Eres — 34mey. However, their 
measurements at frequencies above Eres indicated a strikingly different picture. 
Instead of the gapped spin waves found by other groups, Haydcn et al. re- 
solved four peaks at w > Eres rotated by 45 degrees that dispersed outward 
in energy. Their measured dispersion for spin excitations is consistent with 
that observed earlier in stripe-ordered materials, such as La2-a;Nia;Cu04[79] 
and Lai.875Bao.i25Cu04[61]. Hayden et a/.[64] found little dispersion of the ro- 
tated incommensurate peaks at high energies, 66meV < E < l05meV. This led 
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Figure 18: Smoothed two-dimensional slices through the correlated response 
at 32.5meV and ll.bmeV energy transfers and integrated ±7.b'meV along the 
energy axis in YBa2Cu306.5. At 77.5meV the ring of scattering shows the in- 
tersection of the cone of spin wave dispersion emanating from the (1/2, 1/2) po- 
sition with a constant energy surface. Within statistics the velocity is isotropic. 
The intensity is represented by false color. Reproduced with permission from 
[65]. 



them to conclude that the high-energy spectrum of spin excitations they mea- 
sured was inconsistent with commensurate gapped spin wave spectrum observed 
by other groups. 

High energy spin excitations have also been studied in several materials be- 
longing to the 2-1-4 family. The question of the renormalization of the effective 
exchange coupling J and the spin wave velocity c of the insulator by doped 
carriers was first investigated by Hayden et al.[60] for Lai.86Sro.i4Cu04 in high- 
energy transfer inelastic neutron scattering experiments. They found that the 
spectrum of spin excitations at high energies fits the Heisenberg model well, 
with an effective exchange coupling J that is only mildly reduced from its value 
in the insulator: for this material, Je// = ISOmeV, = 0.15, values that are 
somewhat reduced from the corresponding values in the insulator, J = 153mel^, 

= 0.39. They concluded that quantum fluctuations increase, but that at the 
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Figure 19: Magnetic excitation spectrum and evolution of the form of the mag- 
netic response with energy. Symbols indicate different incident energies, open di- 
amonds: Ei = SOmeV, open squares: Ei = bbmeV, filled squares: Ei = 90meV^ 
filled circles: Ei = 160meV, open circles: Ej = 240meV. " peak-dip- hump" 
structure in the integrated susceptibility x"(^) suggesting that the magnetic 
response has two components. The inverse correlation length k(uj) shows signif- 
icant broadening at high energies. Reproduced with permission from [63]. 
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high energies they study, J and c do not get strongly renormahzed by doped 
holes. However, the spectral weight of their measured high-energy spin excita- 
tion decreases very strongly from its value at x = 0; the peak of that weight 
for X = 0.14 was found[60] to be at w ~ 22meV, significantly below the corre- 
sponding peak for the insulating compound. The details of the spin excitation 
spectrum at low energies, however, could not be resolved. 

More recent high resolution inelastic neutron scattering studies of the stripe- 
ordered compound, Lai.875Bao.i25Cu04[61], and the optimally doped super- 
conductor, Lai.84Sro.i6Cu04[62, 63], reveal the details of the spin excitation 
spectrum at low energies. In their measurements on stripe-ordered sample of 
Lai.875Bao.i25Cu04 Tranquada et al.[61] found two different branches of spin 
excitations that meet at a characteristic frequency — 50 — 55meV, an en- 
ergy spectrum that is similar to that found by other groups, and a sharp fea- 
ture in the spectrum of spin excitations at a; = A2meV. Their results are 
consistent with the behavior seen earlier in other stripe-ordered compounds, 
such as La2-2;Nia;Cu04, where the spectrum of spin excitations was measured 
earlier[80]. In particular, as found by Hayden et a/. [64] and expected on a stripe 
model[61], Tranquada et al. found four incommensurate peaks rotated by 45 
degrees that disperse outward at high energies ui > 50meV. They fit their 
data to a model of charge stripes separated by 2-leg Heisenberg ladders with 
J ~ lOOmeV, a value somewhat reduced from its value J ~ 146mey[57] in the 
insulator. Recent high-resolution neutron scattering studies by Christensenei 
aZ.[62] and VignoUe et aL[63] of Lai.84Sro.i6Cu04 provide further important 
details on the universality of the spin wave spectrum in high-temperature su- 
perconductors. Christensen et al.[62] studied the spin excitation spectrum of 
this material at low frequencies lo < 38meV and concluded that the spectrum 
of incommensurate spin excitations is dispersive inwards and maps onto the 
spectrum of acoustic spin excitations observed in YBa2Cu3 06.85, a material 
with similar planar hole doping. Using this analogy, they concluded that spin 
excitations in Lai.84Sro.i6Cu04 should become commensurate at a frequency 
Eq ~ AlmeV. However, they did not see an associated resonance peak. Vig- 
nolle et al.[63] recently extended these studies to higher frequencies; they found 
the spectrum of spin excitations shown in Fig. 19, with a high-frequency com- 
ponent emerging for lu > iOmeV. Both J and c do get strongly renormahzed 
with doping: the high-energy excitation spectrum they observe corresponds to 
commensurate gapped spin waves with an effective coupling Jeff = SlmeV, a 
value that is dramatically reduced from its value J = lAQmeV in the insula- 
tor. Below u! ~ AQmeV they find the usual[41, 42] attenuated incommensurate 
structure. Their analysis of the unusual double peak structure in the spectral 
weight suggests the spin excitation spectrum separates into two components - 
incommensurate attenuated spin excitations for uj < AOmeV, and commensu- 
rate gapped spin waves for lu > AOmeV. They find that the spectral weight at 
high energies, x"{'^ = l50meV) is roughly 1/3 of the spectral weight observed 
in the parent insulating compound. 

In summary, the energy spectrum of spin excitations as probed by inelas- 
tic neutron scattering is approximately universal, i.e., for different materials. 
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it reflects only the hole doping x. In particular, the positions of incommensu- 
rate peaks at low frequencies are the same in both 2-1-4 and 1-2-3 families of 
materials at similar hole doping levels. The spectrum of spin excitations at all 
frequencies also turns out to be the same at a similar doping level, at least for 
the YBa2Cu306.85 and Lai.84Sro.i6Cu04 pair. The energy spectrum at high 
energies is consistent overall with what is expected for gapped spin waves, and 
measurements on different materials reveal a suppression of the effective ex- 
change constant Jg/ / (x) and the spin wave velocity c{x) with increased doping, 
as well as a suppression of the spectral weight of high-energy spin excitations. 
However, experiments in materials with the hole doping close to 1/8 found a ro- 
tated peak structure at high energies, similar to that expected from static stripe 
ordering. Overall, the energy spectrum of spin excitations in different families 
of materials is consistent with the picture of phase separation and fluctuating 
stripe order that flrst appears below T*. 

2.6 Thermodynamics 

The specific heat measurements and electronic entropy analysis of Loram et 
a/. [16] provide another important constraint on the normal state behavior of 
underdopcd YBa2Cu306+j,. By fitting their data to fermions that are assumed 
to have a "normal state" energy gap of the BCS d-wave form, Loram et a/. [16] 
found that at T = llOif, well above the superconducting temperature, T^, 
their hypothesized gap, A(T, x), increases linearly as the doping is decreased, 
from being ~ zero in a near optimally-doped sample (y = 1) to 200K in an 
underdoped sample with y = 0.7. It is clear from their entropy data that their 
proposed gap size is related to T™-{x) and Jeff{x), that is: 



In a subsequent paper, Loram et al.[17] argued that S/T and the Y Knight shift 
Xs are proportional to each other, with a Wilson ratio for nearly free electrons, 



2.7 Transport measurements 

The peculiar T- linear behavior of the electrical resistivity [81] that has been 
observed near optimal doping in all families of cuprate superconductors, repre- 
sents one of their main unresolved puzzles. It has been explained as representing 
marginal Fermi liquid behavior arising from the proximity of optimally doped 
materials to a quantum critical point[82], with concomitant tu/T scaling[81]. 
However, in the underdoped materials, the electrical resistivity shows strong 
deviations from linearity in T, below a temperature ^ T™. Moreover, scaling 
with a characteristic temperature T™ has been clearly seen in other transport 
measurements, with the best data collapse being that observed in the Hall mea- 
surements of Hwang et al. [20] , who use the following scaling form: 



A(r = llOA') oc T"\x) oc Jeffix). 



(34) 



Xs{T) ~ awj; ■ 



(35) 



i?H = i?^(x) + i?S(x)/(r/r"), 



(36) 
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t=T/r 



Figure 20: The Hall coefficient {Rh) for La2_xSr^Cu04 with > x > 0.34 
plotted rescaled as — R'^{x)]/ R*^ vs t = T/T* . Inset: the parameters 

and i?^ vs Sr composition a;. The temperature T* agrees with T™ obtained 
from the maximum in the bulk spin susceptibility. Reproduced with permission 
from [20] 



where /(T/T™) is a universal function and Rf^{x) and R^{x) arc doping- 
dependent functions (Fig. 20). 

The Hall data for YBa2Cu307_j, from Ito et a/.[19] and Carrington et a/.[21] 
is analyzed in Wuyts et aL[13], who find the following scaling forms (Figs. 21, 
22) for the Hall angle. 9h{T), and the mmrber, n^f(r), lead to excellent data 
collapse: 

cot{9H{T)) = cot{9H{T"^))j^ (37) 

nH(T)=nff(T")^ (38) 

Wuyts et aL[13] extended the T-linear scaling of resistivity near optimal 
doping to the underdoped regime. Their analysis of Ito et a/.'s[19] resistivity 
data for the YBa2Cu307_y family shows that for temperatures T above T™-{x) 
the resistivity is universal and linear (Fig. 23) , and is similar to the optimally 
doped samples: 

p{T)^p{T"^)^ (39) 

The doping dependence of T™ {x) inferred from these scaling forms for the trans- 
port coefficients is consistent with its direct magnetic measurements. 
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Figure 21: Cotangent of the Hall angle cotOn divided by cot 9 h (To) vs square 
of reduced temperature (T/Tq)^ for 15 sets of thin film data together with 
reported data for oxygcn-dcficicnt and Co-doped YBCO single crystals. To{x) 
agrees with the temperature obtained by rcscaling Y Knight shift. Reproduced 
with permission from [13]. 



34 



o YBCO films 
^ YBCO Xtals 
n Co-yBCO xtals 



0.0 I — — ' ' ' ' 1 

05 1.0 1.5 2.0 2 5 

T/T„ 

Figure 22: Hall number nn divided by nniTo) vs reduced temperature T/Tq for 
15 sets of thin film data together with reported data for oxygen-deficient and 
Co-doped YBCO single crystals. To{x) agrees with the temperature obtained 
by rescaling Y Knight shift. Reproduced with permission from [13]. 
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Figure 23: Scaled in-planc resistivity p/p{To) vs scaled temperature T/Tq for 
15 sets of thin film data together with reported data for oxygen-deficient and 
Co-doped YBCO single crystals. To(x) agrees with the temperature obtained 
by rescaling Y Knight shift. Reproduced with permission from [13]. 



Levin and Quader analyzed scaling of the Hall and transport [83, 84] data 
based on the two-band model Eq.(7) with degenerate (^) and nondegenerate 
(77) carriers. They found that the dependence of the cotangent of the Hall 
angle and the deviations from it can be fit to their model if two lifetimes are 
introduced for two different bands, oc and cx T~^. Their two-band 
model predicts three-parameter scaling for the planar resistivity in the form: 

_ A{x)T 
^~ aix)+<i>{W{x)/T)' ^ ' 

where W{x) = r(x — xq) is the energy scale analogous to T*. 

A more recent alternative analysis [22] of the high-temperature Hall data[85, 
86, 87] for the underdoped La2-a;Sr2^Cu04 family has its basis on the idea of 
phase separation. It indicates that the Hall data can be understood if one 
assumes an activated form for the carrier concentration, 

nHaii = no{x) + ni{x) exp {~A{x)/kBT), (41) 

with a doping-dependent no{x) that stays linear in x up to x' 2± 0.12, above which 
strong deviations from linearity are observed. The behavior of A(x) was found 
to be consistent with that measured in photoemission experiments. According 
to Ref. [22], there is crossover behavior when the number of activated carriers 
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becomes approximately equal to the number of doped carriers no(x); this occurs 
at a temperature 

T*{x)^To{x) = -^, (42) 
lux 

that is consistent with characteristic temperatures inferred from other measure- 
ments. At a proposed candidate QCP, x ~ 0.2, the energy gap for activated 
carriers, A(a;), goes to zero[22]. 

2.8 Penetration depth measurements 

One of the earliest universal relations that emerged in the high-Tc field is the 
linear Uemura relation[88] between the superfluid density ps, provided by pen- 
etration depth measurements, and the superconducting temperature, Tc, valid 
in the underdoped regime. This relation provides an important constraint on 
the number of electrons that become superconducting. Recently[89] it has been 
shown (Fig. 24) that a modified Uemura scaling relation, 

PS = 12Q<JdcTc, (43) 

holds for all cupratc superconductors. 

2.9 Angle-resolved photoemission spectroscopy 

Enormous progress in resolution has been made in angle-resolved photoemission 
spectroscopy (ARPES) experiments in the past decade so that ARPES has now 
emerged as one of the best experimental probes of the cuprates[90]. One of the 
most important contributions of ARPES was the detection of the anisotropic 
normal state gap in the fermionic spectrum, first by Marshall et a/. [91], Loeser 
et aL[92], and then by Ding et aL[93] in Bi2Sr2CaCu208+5. 

By examining their data, Marshall et al. concluded that the normal state gap 
has two energy scales. They identified the low energy scale gap 20 — 30TOey) 
by the location of the leading-edge midpoint. This is a clear gap in the fermionic 
spectrum that has a c?-wave-like momentum dependence, with gapless Fermi sur- 
face arcs near the nodal regions. The high energy scale (~ 100 — 200mey) gap 
appears as a broad incoherent feature in the spectrum near the (tt, 0) point, 
where the low-energy spectral weight is strongly suppressed. The low energy 
leading-edge gap first appears below T*[x), a temperature that has the same 
doping dependence as T"^(x) ^ Jeff{x) observed in other measurements. How- 
ever, T* inferred from ARPES mcasurements[90; 28] in Bi2Sr2CaCu208-i-5 turns 
out to be significantly lower than T™, T*{x) ~ T'"(a;)/3 2± J^ff{x)/3. Both 
the measured magnitude and the fc-dcpendence of the normal state leading edge 
gap are very similar to that of the d-wavc gap in the superconducting state; the 
normal state gap smoothly evolves into a superconducting gap below Tc- The 
doping dependence of the normal state d- wave-like gap amplitude A(x) tracks 
that of r*(a;)[90, 28]. 

The ratio of T*{x) and T™(x) varies for different families of materials. De- 
tailed ARPES investigations of the La2_2;Sr3;Cu04[94, 95] and other families of 
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Figure 24: Plot of the superfluid density {ps) versus the product of d.c. con- 
ductivity (odc) and the superconducting transition temperature (Tc) for various 
cuprate superconductors. Odc is measured just above the transition. Within 
error, all of the data fall on the same universal (dashed) line with slope of unity, 
defined by Eq.(43); the dotted hnes are from ps = (120 ± 25)i7£;c?c- Reproduced 
with permission from [89]. 
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cuprate superconductors [90] find the same general form of the spectrum and the 
same hnear doping dependence of the energy gap amphtude A(x) cx T™{x) as 
that found in the Bi2Sr2CaCu208+5 family of materials. Following these and 
other ARPES results, Damascelli et a/. [90] suggest that the value of T*(a;) and 
A(a;) is determined by the maximum value of Tc for a given material, T™"^, 
rather than exchange coupling J: 

A(x) ocT*(x) ocT,""" (44) 

The high energy incoherent feature of the ARPES spectra also has a d-wave 
like dispersion [28]. It is often claimed to be a remnant of the antiferromagnetic 
insulator [90], since it exhibits the same dispersion [28] along the (0, 0)-(7r, 0) and 
(tt, tt) - (tt, 0) directions, which are only equivalent in the reduced antiferromag- 
netic Brillouin zone, and is similar to the ARPES spectra observed in undoped 
antiferromagnetic insulators. 

The fact that the normal state gap closes non-uniformly in momentum 
space with increased temperature or doping has been well established for some 
time[90]. Recently Kanigel et al.[6] investigated the temperature- and doping- 
dependence of the low energy normal state gap in several samples of Bi2212 
and discovered a remarkable scaling relation: the anisotropy of the normal state 
gap and the length of the gaplcss Fermi arcs near the nodal region for these 
different samples depend only on a single parameter, a reduced temperature 
t = T/T*{x) (Fig. 25). Kanigel et al. also found that the length of gapless 
Fermi arcs near the nodes above the superconducting temperature Tc is linear 
in while below Tc the Fermi arcs develop the usual d-wave superconducting 
gap with point nodes. 

Is the superconducting gap that develops at Tc in the nodal Fermi arc region 
the same as the normal state gap, or is it a different gap? Recently Millis[96] ar- 
gued that this important question has not been resolved for years, since different 
experimental studies yield conflicting results. An experimental observation of 
two different gaps in the underdoped regime will confirm that the normal state 
gap in the underdoped cuprates is driven by a different, competing nonsuper- 
conducting order, not pairing fluctuations [97]. Since the superconducting gap in 
the nodal region is extremely small, answering this question in ARPES requires 
a very high resolution study. So far ARPES studies have proved to be conflicting 
as well[96]. For example, a recent leading edge gap ARPES study of Tanaka et 
a/. [98] on three different samples of Bi2212 found that the superconducting gap 
in the nodal arc region had a doping dependence different from that of T*{x), 
and scaled with Tc{x) instead. On the other hand, Valla et a/. [99] reported that 
the energy gap in a non-superconducting sample of Lai.875Bao.l25Cu04, where 
the superconducting order is suppressed by charge ordering, has essentially the 
same simple d-wave form and magnitude as in the superconducting sample at 
higher doping concentration, thus yielding crucial support to the single gap 
scenario. 
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r=55K, r*=500K, f=B.11 

T=100K, r* = 600K, f=0.2 

100 K, r*=26QK, f=D.3e 

10 K, 7" = 240K, f=0.44 

100 K, 7' = 165K, f=0.e 

1S0K, r' = ?40K, f=0.e2S 

ZOOK, 7* = 260K, f=0.76 

140 K, r* = ieOK, f=0.77 

200 K, 7'* = 240K, f=0.B33 

MK, r* = 100K, f=0.9 

1S0K, T*=165K, f=0.91 

2S0K, r'=260K, f=0.96 




Figure 25: (a) Loss of low-energy spectral weight L{(p) = [1 — /(O, ip)/I{A, tp)], 
where /(A, (p) is the symmetrized intensity at the gap energy at the Fermi 
point labeled by angle (p, and /(O, ip) is that at the Fermi energy at the same 
point. Different symbols correspond to different t. The gray area represents the 
straight section of the Fermi surface (b) Variation of the arc length with respect 
to reduced temperature t = T/T*. On the y axis, % is the node and 100 % is 
the antinode. Reproduced with permission from [6]. 
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2.10 STM Experiments 



All STM studies find strong nanoscale inhomogeneity, in which gap magnitudes 
are observed to vary strongly on the nanometer length scale, while the recent 
STM measurements of AUdredge et al. [7] , Gomez et al. [8] , and Boyer et al. [100] 
provide detailed evidence on gap formation in the underdoped 2212 and 2201 
families of cuprates. 

According to Gomez [8], in the overdoped regime locally there is only one 
gap that has a d-wave symmetry. The distribution of the energy gaps is, nev- 
ertheless, very inhomogeneous, indicating the presence of disorder. As the tem- 
perature is raised above Tc, the total spatial area of ungappcd regions increases 
until the whole sample becomes ungapped. The local temperature Tp at which 
the normal state gap first appears varies strongly in space. Nevertheless, the 
ratio of the local gap maximum Ama^. in fc-space to the local ordering temper- 
ature Tp turns out to be universal for different spatial regions and very large [8], 
2^/kBTp = 7.7, much greater than the d-wave BCS limit of 4.3, indicating 
strong coupling. 

In the underdoped regime, Gomez et al. first see the spatially inhomogeneous 
formation of the local gaps below T*, defined as the temperature at which 90 
% of the sample becomes gapped. Below Tc, Gomez et al. observe the forma- 
tion of a smaller kink inside the larger local gap, corresponding to the onset of 
superconductivity. They conclude that in the underdoped regime there are two 
energy scales, with the lower energy scale corresponding to the onset of phase 
coherence. Boyer et a/. [100] adopted a different approach in their analysis of 
their STM measurements on an overdoped sample, a member of Bi2201 family, 
(Bii_j,Pby)2Sr2Cu06+2; (Tc = 15K). Since inhomogeneity is unaffected by the 
onset of superconductivity at Boyer et al. look for the signature of super- 
conductivity, the emergence of the subgap kink at Tc, by removing the effective 
background of the high-temperature STM spectra, which is largely unaffected 
by the onset of superconductivity, from the low-temperature spectra. As a re- 
sult of this subtraction, they find a second, homogeneous superconducting gap, 
Agma;/ = 6.7 ± l.GmeV that forms at Tc — 15K. Their analysis thus strongly 
suggests that the subgap kink at Tc, also seen by Gomez et al., rather than be- 
ing an onset of phase coherence, corresponds to the opening of superconducting 
gap at Tc. 

AUdredge et al. [7] have studied the tunneling density of states in momentum 
space. Their analysis indicates that there is only one gap in the underdoped 
materials. AUdredge et al. find that a large effective anisotropic scattering rate 
r(fc) oc |A(a;| is needed to fit their data in the underdoped regime. Their effec- 
tive scattering rate becomes non-zero at a; < 0.22, and increases approximately 
linearly with decreased doping, similarly to the normal state energy gap. Thus, 
they conclude, they have effectively two kinds of quasiparticles in the under- 
doped regime: quasiparticles near the (tt, 0) antinodcs in momentum space that 
are incoherent and almost localized, and coherent quasiparticles in the nodal 
region. 
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3 A two-fluid analysis of experimental data 



The observations of data collapse reviewed in the previous section provide two 
very important constraints on any theory of the underdoped cuprate supercon- 
ductors. They indicate the presence of two distinct fluids whose composition 
might be expected to vary with doping: a spin liquid containing localized Cu 
spins with a doping-dependent effective interaction and a Fermi liquid whose 
transport properties differ markedly from those of a Landau Fermi liquid. In 
this section we develop a two-fluid framework for analyzing these experiments 
and use it to extract their consequences. We present as well the details of our 
earlier analysis[3] of the magnetic measurements of two fluid behavior. 

3.1 Two-Fluid Description 

Our approach is inspired by the success of the two-fluid phenomenology devel- 
oped for the 1-1-5 family of heavy electron materials[101, 102]. For these and 
other heavy electron materials containing a Kondo lattice of localized /-electrons 
coupled to a conduction band, a two-fluid phenomenological model of the hy- 
bridization of the / electron localized spins with those in the conduction band 
describes very well the emergence of a hybridized non-Landau heavy electron 
Fermi liquid that coexists with unhybridized local f-electrons and conduction 
electrons. The emergence of two components in the bulk spin susceptibility and 
the Knight shift in the cuprates can be understood by writing the total spin of 
the system as a sum of the localized d-clcctron and p-hole spins. 



where are the positions of copper d-electrons, rj the positions of oxygen 
p-holes. Quite generally, the coupling between these spins gives rise to three 
contributions to spin susceptibility [102], 



where Xdd represents the contribution from the localized Cu spins, Xpp represents 
the part of the oxygen p-band that is not hybridized, while Xdp corresponds 
to the magnetic response of the hybridized quasiparticles with a large Fermi 
surface: 




(45) 



X^Xdd + '^Xdp + Xpp, 



(46) 



Xdd 




(47) 



Xdp 




(48) 




(49) 
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The Fermi liquid contribution to spin susceptibility thus arises from Xpd and 
Xpp- Only the position of oxygen holes is important for the above decompo- 
sition; whether or not they form a single hybridized band is irrelevant. The 
three parts are present in an explicit two-band model of Walstedt et al. [36] that 
produces correct expressions for the relaxation rates in terms of Xpp, Xpd, and 
Xdd- It differs from the standard Mila-Rice-Shastry/MMP[32, 40] ionic model 
that assumes a single spin degree of freedom residing on the copper site, so that 
Xpd = Xpp = and the hybridized band description of Millis and Monien[103], 
in which although Xppi Xpdi 

and Xdd are taken to be non-zero, all three com- 
ponents are assumed to track the spin susceptibility of a single hybridized pd 
band. 

The scaling results of the previous section tell us that Xdd must maintain its 
local character throughout much of the phase diagram, for how else could the 
bulk susceptibility and the low frequency magnetic response map onto the 2D 
Heisenberg model? To describe this we therefore write: 

Xdd = f{x)xsL{q,i^) (50) 

where f{x) is the fraction of the dd (and the total) response function that retains 
its local spin character, while the remaining portion of the spin-spin response 
function describes the response of a (non-Landau) Fermi liquid with the large 
Fermi surface that results from hybridization: 

2Xdp + Xpp = [1 - f{x)]xFL{q, t^)- (51) 
We thus arrive at the two fluid description, 

x(g, = f{x)xsL{q, + f{x)]xFL{q, w), (52) 

that we will use in our subsequent analysis to extract from experiment the 
doping dependence of both xsl and f{x). 

The spin liquid contribution completely determines the low frequency dy- 
namic magnetic spin susceptibility measured by NMR on copper nuclei. It 
corresponds to the scaled spin response of the 2D Heisenberg model for local- 
ized copper spins, which in the low frequency limit can be written in a general 
phenomenological form proposed by Millis, Monien, and Pines[40], modified to 
include the possibility of propagating spin wave excitations [104]: 

Y5i(<7,tj) = 5—, (53) 

i + eiq-QY -I— - ^ 

where Ag^ = c/^ is the gap in spin excitation spectrum. We note that in general 
damping in the Heisenberg model is caused by the spin wave scattering, so that 
the relaxational frequency, 0^5^, must be, in general, frequency-dependent in the 
underdoped cuprates. We further note that if one assumes Heisenberg model 
scaling, then, from dimensional arguments, a{x) = XQ^~^ is proportional to 
l/T'^[x). 
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The Walstedt et al. generalization of the Mila-Rice-Shastry hyperfine Hamil- 
tonian for Cu site can be written as: 



Here i corresponds to the nearest neighbor Cu sites, p = a,b,c while D is a new 
transferred hyperfine constant for the spins sitting on oxygen. The sum index j 
goes over the oxygen sites neighboring the copper site. Similarly, one can write 
the hyperfine Hamiltonian for the oxygen nuclei: 
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hyp 



^|"/,c5:S,^(r,)+"/pi?pS^| (55) 



Here C is the transferred hyperfine coupling constant for localized Cu spins and 
Ep is the contact hyperfine interaction for the oxygen nucleus; the sum index i 
goes over the NN copper sites. Making use of 

fiB{SP{r)) = ixpd + Xpp)Be.t (56) 
fisiS'^ir)) = {xpd + Xdd)Be^t, (57) 

where Bf-^t is the strength of the external magnetic field and /is is the electron 
magnetic moment, we get, from the above Hamiltonians: 

63A'p = ^^IU,p + ^^^^ ixdd + Xpd) + g3 2 iXpd + Xpp) (58) 

2C E 

1% = ^'Ko^p + - (^^rf + ^p^) + _ ^(^p^ + ^^p). (59) 

Here 7„ are the corresponding nuclear gyromagnetic ratios. On introducing the 
dynamic spin susceptibilities, 

X^.AQ,^)^{{S'p{q)S£{-q)))^. (60) 

where fi,v = p, d, we obtain the Walstedt et al. expression for relaxation rates: 



Jmx^uiq,^^) 



li- E "^...(9)^^^^^^^, (61) 



,2 >^2 



where 7 is the direction of magnetic field, a is the nuclear site. For Xdd one gets 
the usual MMP form factors: 



^^^||,dd(9) = iA± + 2B{cos{q,a) + coa{qya))\ (62) 

""^F^^q) ^^e^F^^^Miq) + ^ {A\\ + 2B (cos q^a + cos qya)f , (63) 
^'F^Mq) = 2C2(1 + cosiq^a)). (64) 



44 



The form factors for the mixed p-d contribution are: 



\\,pd\<i) - ^\\.dp\ 

2i?(cos {qxa/2) + cos {qya/2)){A±^ + 2B{cosqxa + cosqya)), (65) 



= D(cos (g^a/2) + cos {qya/2)){A± + + 4:B{cosqxa + cosg^a)), (66) 

''^^11,^^(9) J^||,dp(9) = 2C^i cos (g.a/2), (67) 

i^Fi,pd(q) ^^±,dp(g) = C(£;i +S||)cos(g,a/2). (68) 
Finally, for the p-p part of spin susceptibility: 

^^F^,pp(q) = 4i:>2(cos ((?:ra/2) + cos {qya/2) f (69) 

''i^ll, pp(g) = i^i (70) 

i^^^^,pp(9) = ii?i + i4 (71) 

Apart from ^'^Fpp{q), the form factors for the p-d and p-p contributions 
vanish at Q = (tt/o, 7r/a): 

^'FpdiQ) ^""^ FpdiQ) FppiQ) = 0. (72) 



In what follows we focus on the behavior of planar dynamic and static magnetic 
response and the thermodynamic behavior, as measured by the inelastic neutron 
scattering, NMR, heat capacity and bulk spin susceptibility. 



3.2 Scaling for bulk spin susceptibility and Knight shift 

We now apply our two-fluid description to an analysis of the bulk spin suscep- 
tibility and Knight shift measurements in La2-xSr2;Cu04 and YBa2Cu3 064-x 
families of cupratc superconductors. For the static susceptibility, the two com- 
ponents in Eq.(52) can be written in the following form: 

xiT) = f{x)xsL{Tir^{x)) + (1 - f{x))xFL + Xo, (73) 

where 

Xo = XVV + Xcore- (74) 

Here xvv and Xcore are temperature- and doping-independcnt Van Vleck and 
core contributions. Xo is usually taken to be diamagnetic, due to a large core 
contribution[2]. As noted earlier, xsl follows very well the calculated[5] bulk 
spin susceptibility for the Heisenberg model with a doping-dcpendent exchange 
constant Je//(x) T™(x). Thus, the spin liquid contribution to static magnetic 
response can be written as 

f{x)xsL{T/T"^{x)) = x'"x(T/r™(x)), (75) 
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Figure 26: Comparison of the scaling curve for bulk spin susceptibility data in 
underdoped metallic La2-2:Sr2:Cu04 from Nakano et al.[ll] to the Heisenberg 
model numerical calculations of Makivic and Ding[5]. The maximum for spin 
susceptibility is reached at T™ ~ 0.93 J. Deviations from the Heisenberg model 
results are observed for T < T* ~ r'"/3. The inset shows the numerical 
results for the correlation length[5], that demonstrate that ^ ~ 1 at temperature 
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where x™ is the maximum value of the spin liquid susceptibility and x(T/T'"(x)) 
is a universal function of T/T™. A comparison of the universal data collapse 
curve for the bulk spin susceptibility in La2-a;Sr2;Cu04 family found by Nakano 
et a/.[ll] to the Heiscnberg model calculations [5] is shown in Fig. 26. We see 
that the universal function of Nakano ei al. deviates significantly from the 2D 
Heisenberg model results at low temperatures, T <T* ~ r™/3. 

As noted earlier, the scaling analysis of Nakano et al.[l\\ differs from that of 
Johnston[2], because it includes at lower doping levels a Curie term C/T, and a 
linear term B{T — To) in the static spin response. While the inclusion of these 
terms may be justified, it introduces a number of new parameters that are, in 
a sense, unnecessary since the universal curve in the insulating compound must 
agree with the 2D Heisenberg model results. We therefore omit these terms in 
our analysis, and we assume, in agreement with recent photoemission studies [6], 
that the Fermi liquid contribution, xfl(T, x), is temperature- independent above 
T* , but could, for the 1-2-3 materials become temperature-dependent below the 
temperature, T*, at which a gap starts to form in the quasiparticle energy 
spectrum, leading to the formation of the Fermi arcs. 

To determine the doping dependence of the parameters in Eq.(73), we ana- 
lyzed the bulk spin susceptibility data[ll] for metallic underdoped La2-xSr3,Cu04 
With a constant contribution Xo+X-Fl(^) subtracted, we find, in agreement with 
the earlier results[2, 12, 11], the data collapse to the 2D Heisenberg model curve, 
shown in Fig. 27, where we see that the data collapse continues below the tem- 
perature ~ T^/S at which one no longer finds agreement with the Heisenberg 
model. 

A similar scaling analysis for the ^^Cu Knight shift data in 1-2-3 materials is 
shown in Fig. 28. The doping dependence of T"^{x) ~ 0.93Je//(a;) that follows 
from the analysis of the bulk susceptibility and Knight shift data, is shown 
in Fig. 29 for both La2-xSra;Cu04 and YBa2Cu306+x families. T™{x) falls 
linearly with increased doping x for both materials. For La2-xSra;Cu04, and 
doping levels less than 0.18, 

r™(a;) = 1218ii:(l - 4.45a;). (76) 

Our analysis of the bulk susceptibility data in La2_2^Sr^Cu04 also enables us to 
extract the doping dependence for the Fermi liquid and spin liquid components. 
The results for x™{^) ^-nd Xfl{x) given by Eqs.(75),(73) are shown in Fig. 30. 

The doping dependence for the fractional occupation, f{x), of the spin liq- 
uid state can be determined from T"^{x) and x'"(x). The dimension of x is 
1/ energy. Thus, since the universal curve for x > T* coincides with the nu- 
merical calculation for the 2D Heisenberg antiferromagnet with T™ ~ 0.93Je//, 
we can write: 

X5l(t/t™(x)) a x(r/r™(x)), (77) 

where x(T/T"^{x)) is a universal function for the Heisenberg antiferromagnet. 
Using Eqs(77) and (75), the doping dependence of f{x) is then determined by 
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Figure 27: Comparison of the bulk spin susceptibility data in underdoped metal- 
lic La2-a;Sra;Cu04 from Nakano et a/. [11] to the Hcisenbcrg model numerical 
calculations of Makivic and Ding [5]. The maximum for spin susceptibility is 
reached at T™ ~ 0.93J. While the data collapse remains reasonably good even 
at low temperatures, large deviations from 2D Heisenberg model prediction are 
observed for T < T"'/3. 
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Figure 28: Scaling for the ^^Cu Knight shift[37, 33, 105, 106, 107, 108] in 
YBa2Cu306+a:, Compared to the scahng function obtained by Nakano et aZ.[ll] 
for the bulk spin susceptibility in La2-a:Sra;Cu04. T™ is the temperature at 
which the Knight shift has a maximum. 



the product x'"(2:)r™(x): 

X™(-^ = 0)r™(a: = 0) ^ ' 

From our data analysis, we find that for the doping levels available for our 
analysis, f[x) also decreases linearly with x, 

/(x) = 1 - 4.977.T, (79) 

as shown in Fig. 31, so that T"^{x) and f{x) roughly track each other. Fig. 30 
also shows the doping dependence for the Fermi liquid component, xfl{x) = 
(1 — f{x))xFL + Xo- We find that the Fermi liquid component increases linearly 
with doping in the metallic regime, in agreement with Eq.(73). The offset for 
this linear dependence is due to the presence of xo: which is diamagnetic [2]. 
Since f{x) also decreases linearly with x, we conclude that xfl in Eq.(73) has 
only mild, if any, doping-dcpendence. 



3.3 Knight shift data 

We now examine the Knight shift data in more detail. As noted in Section 
2, most early Knight shift experiments on different nuclei found the same 
temperature-dependent contribution to spin susceptibility [33, 34, 109], although 
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Figure 29: (a) The temperature T"^{x) that follows from our analysis of the 
Nakano et a/. [11] bulk spin susceptibility data for La2_a;Sra;Cu04, and (b) the 
Knight shift data for YBa2Cu306+2; 
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Figure 30: x™(a;) and XFLix), from our analysis of the bulk susceptibility mea- 
surements in La2-2;Sr2:Cu04. Here xfl{x) here also includes a dianiagnetic 
constant, xo, hence the non-zero offset. 




0.00 0.05 0.10 0.15 0.20 

doping, x 

Figure 31: The spin liquid fraction f{x), from our analysis of the bulk suscep- 
tibility measurements in La2-2;Srj;Cu04[ll]. 
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some deviations from this universality (within error bars) can be seen on the 
universal plot below the superconducting temperature Tc [33, 34]. The linear 
K — X and K ~ K plots above Tc, observed in many early Knight shift experi- 
ments, then provide a measurement of the hyperfine couplings to single copper 
spin. These measurements had been regarded as a key proof of an effective one- 
component model. However, the Cu Knight shift measurements on the 1-2-4 
material by Suter et aL[38] and the recent apical oxygen Knight shift measure- 
ments in Lai.85Sro.i5Cu04 by Haase et aL[39], tell a different story, and we 
consider these now. 

The Knight shift measurements of Suter et al. were done for an external 
magnetic field in the c direction, where, because of the well-known accidental 
cancellation [40] of the Xdd form factor, A|| + AB, in Eq.(58), within the single- 
component approximation, the Knight shift for fields along the c-axis would be 
of purely orbital origin. Within the two-fluid description, we see from Eq.(58) 
that one can, in addition, probe Fermi liquid behavior through an isotropic 
transferred hyperfine interaction D that arises from the hybridization of planar 
oxygen p-orbitals and copper s-orbitals. So to the extent that the Fermi liquid 
susceptibility becomes temperature dependent as a result of a quasiparticle gap 
opening up at T*,one would expect to see this in As may be seen in 

Fig. 32, where we have extracted the behavior of Xpp + Xpd as a function 
of T from their Knight shift measurements, this is just what was found by 
Suter et a/. [38]. The anomalous Fermi liquid behavior begins at a temperature, 
T* ~ 180K ~ 0.4T'" for YBa2Cu308, an onset temperature that is consistent 
with the onset of Fermi arc behavior at a comparable value of T* seen in the 
ARPES measurements on underdoped 2212 materials. 

In their recent re-analysis of the planar ^^Cu and ^^O Knight shift experi- 
ments on Lai.85Sro.i5Cu04 that was accompanied by new apical oxygen Knight 
shift experimental results, Haase et aZ.[39] have found convincing experimental 
evidence for the existence of both a Fermi liquid and a spin liquid component 
in this material. They show that while all nuclei see the same temperature- 
dependent spin liquid component, the different nuclei see as well a component 
that is independent of temperature above the superconducting transition tem- 
perature, Tc, but displays the expected d-wave signature below Tc, and hence 
must be associated with a Fermi liquid component, rather than representing an 
orbital or Van Vleck core term. 

Quite importantly, they combine their results in the superconducting state 
(that include correcting for any Meissner shielding effects) and the normal state 
for magnetic fields in and perpendicular to the CuO planes with Equations 
(58) and (59), supplemented by an equivalent expression for the apical oxygen 
site, to determine all the relevant hyperfine couplings. They show these can 
be used to determine the strength and temperature dependence in the normal 
and superconducting states of all three components, Xdd^ Xdpi and Xpp- 

We 

refer the interested reader to their paper for the details of their findings, which 
place important constraints on any description of the hybridization of the Cu d 
electrons with the oxygen p band. 

While further Knight shift measurements under conditions favorable to ob- 
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Figure 32: Normalized temperature dependence of Xpd as a function of T/T"^ , 
inferred from ^^K^ measurements [38] in YBa2Cu308. 



serving the Fermi liquid component are highly desirable, we believe it is rea- 
sonable to conclude that the two-fluid analysis of Knight shift in heavy electron 
materials[102] extends naturally to high-Tc superconductors, with 

Xdd = f{x)xsL, (80) 

and 

Xpp + 2Xpd = {l- fix))xFL. (81) 

Eq. (81) can then be combined with the results of Haase et ai, and our previ- 
ously determined value of f{x), to determine the doping- independent quantity 
Xfl for the 2-1-4 materials. The results of this analysis can be checked by 
extending the Haase et al. analysis to other doping levels. 

3.4 63 Cu NMR Relaxation Rates 

In analyzing the NMR relaxation rates, we shall focus on the ^^Cu spin-lattice 
relaxation rate, ^"^Ti, and the spin-echo-decay time, ^^T2G- As discussed in the 
beginning of the section, there are different components of the spin response. 
However, due to the large peak of the spin-liquid part at Q = (7r,7r), and the 
corresponding non-vanishing form factors for ^^Cu nuclei, the copper relaxation 
rates primarily probe the low-frequency properties of the spin-liquid contribu- 
tion that arises from Xdd = f{x)xsL part of the spin susceptibility. In our 
modification of MMP[40] both Xpd and Xpp correspond to the contribution of 



53 



fermions with a large hybridized Fermi surface, or pieces of such a Fermi sur- 
face, in agreement with the Knight shift measurements of Haase et al.[39], and 
therefore have a neghgible effect on copper relaxation rates. An alternative pos- 
sibility, that we mention only briefly, is that the two components seen in the 
bulk measurements are a result of microscopic phase separation and formation 
of dynamic stripes. Surprisingly, while the transition points that we obtain from 
the two-component analysis of bulk measurements coincide with the well-known 
region of phase separation in La2CuO4+5[110], the two-component nature is ev- 
ident in bulk measurements at temperatures well above the phase separation 
temperature region. 

Where these can both be carried out, measurements of ^^Cu spin-lattice 
relaxation and spin-echo decay rates indicate the presence of two different dy- 
namic scaling regimes that lie above and below T^IQ]. Above T™ one is in a 
mean-field, z = 2 scaling regime, in which the relaxational frequency, ujsf varies 
as 1/^^; below T™ it varies as 1/^ because the spin hquid has entered the z = 1 
dynamic scaling regime[9, 18] expected for the quantum critical (QC) regime 
of a 2D antiferromagnet. As discussed in Ref.[9], one can describe the spin dy- 
namics using the quantum non-linear sigma model, or spin wave theory [47, 48]. 
The resulting QC scaling theory[48] for the spin liquid without long-range order 
gives a linear dependence of the correlation length on temperature, 

^ 1.04- +aix), (82) 



e(T,x) • c 

where c (x J ^ T™ is the spin wave velocity, and the offset a{x) > goes to 
zero at Xc, a quantum critical point for the spin liquid that marks the onset of 
long-range order. A similar linear dependence on T can be expected also for 
63riT oc cosf[9, 18], 

e^TiT = A{x) + ^T™^ , (83) 

with an offset A{x) that measures the distance from the proposed quantum 
critical point, while k is a universal coefficient. 

Quantum critical theory applies in the universal regime, ^ ^ 1, thus, well 
below T"'. Scaling of the form seen in Eq.(83) is self-evident from the low- 
temperature ^^Ti data on La2-a:Sra;Cu04 material of Ohsugi et aL[49](Fig.7). 
When they plot the product ^^TiT vs temperature in the metallic undcrdoped 
regime, they find a set of parallel lines for their different doping levels that 
extend from T ~ 300K down to a low-temperature upturn in ^^TiT that is near 
Tc- We plot in Figs. 33, 34 the scaling behavior for ^'^Ti in the 2-1-4 and the 
1-2-3 materials. The x dependence of the off-set, A{x), of the low-tcmpcrature 
linear behavior in ^^TiT points to a critical point in the spin liquid at x = 
0.05, separating long range order from short range order, a result that was also 
suggested earlier in an analysis of the high-temperature ^^Ti data[18]. The 
long range magnetic order ai x < 0.05 can be a 2D antiferromagnet or a spin 
glass [3]. 

While the details of the microscopic theory at low doping may vary, the ex- 
perimental data imply the presence of a linear spin wave-type excitation spec- 



54 





0.08 -| 




0.07- 




0.06- 
0.05- 


< 


0.04- 




0.03- 


1- 

« 


0.02- 



0.01 



T = A(x) + T"F(T/T") 
F(T/T'")=C in'" 



• 


x=0.075 


■ 


x=0.1 




x=0.13 


♦ 


x=0.15 


T 


x=0.2 



.t 




^ 20 
U 

S 15 




QC point 



— I — 
50 



100 150 200 250 

T,K 



doping, X 
300 



Figure 33: Scaling for the NMR relaxation rate Ti in La2_xCu04[49] The offset 
A{x) shown in the inset depends linearly on x, and points to a QCP in the spin 
liquid. 
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Figure 34: Scaling for the NMR relaxation rate Ti in YBa2Cu306+2;. The offset 
A{x) shown in the inset depends linearly on x, and points to a QCP in the spin 
liquid. 
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trum in the spin liquid with a temperature-dependent gap, Asl{x,T), that is 
controlled by its distance from the quantum critical point at x = 0.05. The spin 
wave excitation spectrum has the following form: 



6(g) = ^A|^ + c2(g-Q)2, (84) 

Here c is the spin wave velocity. In the z = 1 QC regime the spin gap Asl 
tracks ujsf = c'/^. For x > 0.05 the gap saturates at low temperatures to 
Ao = c/^o{x) at the crossover to the QD gapped spin hquid regime. We note 
that the origin of the saturation of the correlation length is largely irrelevant for 
data analysis - an energy gap can appear in the 2D quantum spin liquid[47, 48], 
or in ID stripes of spin liquid[61] due to dimensional crossover. In the latter 
case the correlation length will saturate at the size of the domain or dynamic 
stripe order, ^(T = 0) oc i. Therefore, experiment docs not distinguish between 
these possibilities. 

The QC-QD crossover could explain the sharp upturn in ^^TiT at low 
temperatures observed in NMR experiments on 1-2-3 materials at tempera- 
tures well above Tc that is not found in experiments on the 2-1-4 materials. 
Since c{x) oc T'"(x), and becomes small as one approaches x = 0.22, while 
1/(^q(x) increases linearly from x = 0.05, the zero-temperature spin liquid gap 
Ao = Asl{x,0) will have a bell-shaped form similar to Tc{x). The gapping of 
the spin liquid leads to an exponential decrease of damping with temperature 
in the spin liquid. 

The detailed properties of the spin liquid can be extracted directly from the 
NMR ^^Ti relaxation rates measurements using the ansatz for the correlation 
length ^(T™) = 1, or, where these exist, from both ^^Ti and ^^T2g measure- 
ments, as it was shown some time ago by the authors[9]. In the z = I QC 
regime, 

''T,T oc (85) 
a{x)f{x) 

^ oc a(x)/(x)e, (86) 



where all parameters are taken from Eq.(53); these equations differ from those 
in Ref. [9] by the presence of an extra factor f{x) in Eq.(53). A combination of 
these two measurements [9] gives a handle on the spin wave velocity c: 

^5^7— OC wsfC = c' c:i 0.55c (87) 

The coefficients of proportionality are given by the hyperfine Hamiltonian. As 
easily seen from the above equations, 

63ri(r™)r™oc— (88) 

f{x)a{x) 

The high-temperature scaling for ^^Ti, ^•^Ti(T) = const, that was seen by 
Imai et a/.[18] in La2-i^Sr^Cu04, but is not present in the 1-2-3 materials. 
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Figure 35: c'{x) from NMR experiments in YBa2Cu306+: 



at first sight seems to imply x- independent exchange integral, J, or c', which 
contradicts the bulk susceptibility data for this material. However, a closer look 
shows there is no contradiction. As we have seen in Section 3.1, f{x) and T"^{x) 
both depend linearly on x. Thus, on general grounds, one can write: 

c'{x) oc c{x) oc T™(x) oc a~'^{x) oc f{x). (89) 

The Imai et al.[18] scaling then follows from Eq.(88), since one can approxi- 
mately write: 

a{x)f{x) ~ a(0) ~ ll.3states/eV (90) 

In general, the empirical scaling of Eq.(89) is somewhat accidental, since the 
low-energy properties of the spin liquid can be renormalized by interaction with 
fermionic quasiparticle excitations. This renormalization appears to be small 
for the 2-1-4 materials but appears to be important for the 1-2-3 materials, 
where the Ohsugi et a/. [49] results and the Imai et a/. [18] high-temperature 
scaling of ^'^Ti are not observed. f{x), T™{x), and c'{x) can in principle be 
independently determined from NMR experiments. Since ^^T2g has only been 
measured for the YBa2Cu306+x family, we plot the doping dependence of c' 
for these materials in Fig. 35, assuming commensurate local spin response. It 
is evident that c'{x) decreases approximately linearly with doping. However, a 
comparison of c'{x) with T™'{x) shows that the two quantities do not track each 
other, violating the empirical scaling of Eq.(89). While the doping dependence 
for both is linear, it is evident from Fig. 36 that they are not proportional to one 
another. The violation of the Imai et a/. [18] scaling relation by this family of 
materials is therefore to be expected. 
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The doping dependence of the parameter c' {x)/ f{x)a{x), that follows from 
NMR ^^Ti measurements only, is shown in Figs 37,38. We see that the empirical 
scaling relation Eq.(90) holds reasonably well for the La2_a;Sra;Cu04 family, so 
that the doping dependence of this parameter also yields the doping dependence 
of c(x). 

3.5 Thermodynamics 

We now turn to an analysis of thermodynamic measurements of the electronic 
entropy[16]. These are quite difficult to carry out, since the phonon contribu- 
tion accounts for 98-99 % of the measured value of the entropy, and must be 
subtracted to a high degree of accuracy to get the electronic contribution. The 
subtraction procedure is a complicated process, involving subtracting data for a 
reference material, such as YBa2Cu306[16], and making additional corrections 
for the modification of the phonon DOS with doping. The reference material 
is usually an insulator, either the parent compound, or a compound in which 
the Fermi liquid is suppressed with impurities such as Zn. The difficulty with 
such a subtraction procedure for analyzing the data is that the spin liquid con- 
tribution exists even when the Fermi liquid is suppressed. Thus, a consistent 
data analysis requires that the reference material be known. For this reason 
we restrict our attention to the data of Ref.[16], where the reference material 
is YBa2Cu306. The spin liquid contribution of the parent compound is sub- 
tracted from the electronic entropy in this analysis. Nevertheless, at moderate 
temperatures T < 300-ftr it remains comparatively small. 
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Figure 37: c'{x)/ f{x)a{x) from NMR experiments in La2-xSra;Cu04 
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Figure 39: Scaling for the spin liquid part of the entropy Ssl{T /T™-) in 
YBa2Cu306+a;. The entropy data is from Ref.[16] in the intermediate doping 
regime, where the spin liquid part of the entropy is not subtracted. T™ (x) and 
/(x) obtained from this analysis are consistent with our analysis of the NMR 
Knight shift and the bulk spin susceptibility. The solid line is the integrated 
numerical result for the heat capacity of the 2D Hcisenbcrg model taken from 
Ref. [5]. 



Our analysis of the data of Ref. [16] shows that the entropy displays scaling 
and two component behavior in the form: 

S = /(x)55L(r/r") + (1 - /(x))7T, (91) 

with a scaling function Ssl{T/T™') that corresponds to the theoretically calcu- 
lated entropy for the Heisenberg model [5]. We have integrated the numerical 
results for the heat capacity of the Heisenberg model from Ref. [5] to get the 
expected entropy. We show the results of this analysis and the data collapse 
on the Heisenberg scaling curve in Fig. 39. We find that at small doping scaling 
becomes progressively worse, as might be expected because of the subtraction 
procedure followed. 

To understand our result, we note that on general grounds, one can write 
the entropy in the form: 

S = .f{x)SsL + (1 - fix))SFL + Sconf.g, (92) 

where Scon fig is the configurational entropy of the fluctuating stripe or domain 
order discussed in the next subsection. The configurational part, however, is 
not seen in the data factorization Eq.(91) We believe, the reason for this is that 
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the configurational contribution of the fluctuating stripes could also be linear 
in temperature[lll]. Indeed, the domain walls can be treated as solitons, with 
the Hamiltonian 

1 f , f du 



where y is the direction along the stripe, while u is the displacement of the 
wall perpendicular to its equilibrium direction; eq is the linear tension of the 
domain wall (soliton). This model maps on a ID fcrmion modcl[lll, 112], and 
the resulting configurational entropy has a linear form, 

Sconfig = -—j2 i (94) 

where I is the mean distance between solitons. Such a linear contribution to the 
entropy is indistinguishable from the Fermi liquid contribution. 



3.6 Inelastic neutron scattering 

Our analysis also casts light on the issue of intrinsic spatial inhomogcncity[53] 
in the pseudogap state[79, 61]. The issue is how to reconcile the commensu- 
rate spin susceptibility peaks required by the spin liquid and the fit to the 
NMR data[36], with the incommensurate peaks seen in inelastic neutron scat- 
tering experiments [26]. The oxygen relaxation rate is known to be very small 
compared to copper, with a very different modified Korringa-type temperature 
dependence. An attempt to explain it within the one-component model with 
incommensurate spin liquid peaks[36, 45] in Lai.85Sro.i5Cu04 leads to an oxy- 
gen relaxation rate that is a factor of three too large compared with actual 
experimental data, and has the wrong temperature dependence. 

Similar inconsistencies occur with relaxation rates on other nuclei, such as O 
or Y in YBa2Cu306+a;. The contradiction cannot be resolved if one assumes two 
different components sitting on oxygen and copper. Indeed, since the Knight 
shift on different nuclei has an identical temperature dependence[34, 39], it 
requires the existence of a transferred hyperfine coupling from the copper [32]. 
Such a transferred hyperflne interaction assumed by MMP will, however, induce 
a "leakage" spin liquid contribution from the spin liquid in the two-component 
picture for any incommensurate spin liquid response. 

The resolution was suggested some time ago by Slichter[44], who pointed 
out that if, as a result of Coulomb interaction between the holes, regions that 
were hole rich formed domain walls between regions in which one had nearly 
antiferromagnctic commensurate behavior, the two kinds of experiments were 
compatible. This issue has been independently confirmed and clarified by Tran- 
quadaei aL[79, 61], who found that the spin wave dispersion seen in cuprate 
superconductors by inelastic neutron scattering is consistent with commensu- 
rate spin response of a spin ladder, while the incommensurate structure appears 
as a result of the fluctuating stripe order. 

Such a scenario implies finite-size scaling, and saturation of the antiferro- 
magnctic correlation length at low temperatures at a value of the order of the 
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size of the domain. Estimating the size of the domain from the discommensu- 
ration seen in neutron scattering experiments [26] in Lai.86Sro.i4Cu04, Sq = ttS, 
S = 0.245, we can write: 

SqL = 271 = tt6L (95) 

Thus, 

L/a=^ = 8.16, (96) 



where a is the lattice spacing. Acppli et aL[26] find the saturated correlation 
length in La2_a:Sra;Cu04 at a doping level of x = 0.14 is 

Co = 29.41 = 7.7a, (97) 

a value that is very close to our estimate of L. We conclude that both magnetic 
probes imply the presence of a fluctuating striped or other domain order. Apart 
from the issue of incommensurate peaks seen in INS experiments and its possible 
resolution in terms of a fluctuating domain order, the spin wave spectrum and 
the correlation length observed in inelastic neutron scattering experiments is 
completely consistent with our proposed two-fluid scenario, as we show below. 

We return now to the magnetic scaling seen in the inelastic neutron scatter- 
ing experiments which, like the ^^Cu NMR relaxation rates, probe spin liquid 
behavior. We begin with the experiments of Sternlicb et aL[25] that imply the 
presence of a true uj/T scaling for the local spin susceptibility in the underdoped 
material YBa2Cu306.6- The data collapse observed experimentally can be fit 
to the following simple scaling form: 

XU^,T) = A{x)^ , (98) 

with deviations from scaling that occur at progressively lower temperatures, as 
the frequency lu increases. 

This scaling confirms the high-temperature scaling seen in NMR[18]. To 
clarify this point, we plot both experimental results in the same units in Fig. 
40. The quantity measured in neutron scattering experiments analogous to 
^^TiT oc OJSF is a function of the local spin response, u/x'ii^^T), typically 
given in counts/minute. We therefore fix the absolute value of that response to 
be consistent with high-temperature NMR. Both experiments reveal linear in 
temperature cjsf, a signature of magnetic quantum critical behavior. At finite 
frequencies this regime extends to lower temperatures, as one may expect on 
general grounds. 

The presence of z = 1 quantum critical behavior is verified further by the 
INS measurements of the temperature dependence of the correlation length seen 
in Lai.86Sro.i4Cu04[26]. Aeppli et aL[26] show that 

^,=^^+^[{kBTf + {n^)\ (99) 

where = 29. 4A, V ~ A7meV. For the QC scaling in the cr-model[47, 48], 

I = 1.04fcBT. (100) 
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Figure 40: Local spin response seen in inelastic neutron scattering in 
YBa2Cu306.6, w/xi(w,T) , scaled to NMR ^^TiT in La2-:ESr:rCu04 and 
YBa2Cu306+2:. 



On the other hand, the Aeppli et al.[26] result gives, if ^ is measured in lattice 
constants: 

(101) 



Thus, for X = 0.14 we find 



c= 1.04r = 49mey, (102) 



a value that is much less than that found for the insulator value Cma = 22AmeV, 
and is qualitatively consistent with the behavior of the effective J inferred from 
the bulk spin susceptibility measurements. 

Another comparison of theory with experiment on the spin liquid phase 
diagram comes from the neutron scattering measurements of the magnetic cor- 
relation length. For the Heisenberg spin liquid, the correlation length ^ = 1 at 
T ~ T™. On using the results for ^ obtained by Aeppli et ai, we find: 

= 1) = T"{x = 0.14) = 540K ± lOOK. (103) 

In good agreement with the result, T™ = 460K, a,t x = 0.14 doping obtained 
using Eq.(76). 

A further important point that follows from the results of Aeppli et al.[26] 
concerns the breakdown of scaling, which should occur at a universal scale Ta, 

such that 

keTa = Vi^^)^ + {ksTU^W. (104) 
where is the temperature of the breakdown oiuo/T scaling seen in the neutron 
scattering experiments at a fixed frequency u. NMR experiments just measure 
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Figure 41: Temperature-frequency analysis of scaling failure, data from Stern- 
lieb et al. and Takigawa (6.63 NMR) 



Ta directly. We analyze the data of Sternlieb et al. [25] using this framework to 
see if it is consistent with the results shown in Fig. 41. According to Aeppli et 
al., the breakdown of scaling frequency/temperature can be fit to the form: 

= El - ATI (105) 

From this analysis we find, for YBaCuOe.e: 

Ebr = 8.03mey, A = -0.00256{meV^ / K^), (106) 

or 

S.OSmeT/ = 93K, (107) 

159K, (108) 

0.58. (109) 

Overall, this experimental picture is self-consistent, although the experimental 
finding of the simple form for the breakdown Eq. (105) is rather surprising. 

We conclude this section with a brief discussion of phase diagram for the 
spin liquid that is consistent with the experimental constraints discussed above. 
Our proposed phase diagram[3] is shown in Fig. 42. 

The phase transition or a crossover at a doping level ^ 0.22 is that from 
the conducting Fermi liquid on the right-hand side to pseudogap matter on the 
left-hand side, in which a portion of the quasiparticle Fermi surface has lost its 



UJbr[T=Q) = 

Tbriu; = 0) = T"'/3 = 

Thr 
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Figure 42: (a) The spin liquid phase diagram for La2_2;Sr2;Cu04. T™(x) corre- 
sponds to a crossover at ^(T, x) = 1 to the strongly correlated phase in the spin 
liquid. T*{x) ^ T™(a;)/3 corresponds to the appearance of fluctuating stripe 
order Tgap marks the spin liquid crossover from the quantum critical (QC) to 
quantum disordered (QD), or gapped spin liquid, regime in the stripe phase. 
Tqc-rc is the spin liquid crossover from the QC to renormalized classical (RC) 
regime, (b) The spin liquid phase diagram for YBa2Cu306+i: displays similar 
crossovers. 
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low-frequency spectral weight, becoming localized in an insulating spin liquid 
that coexists with the remaining Fermi liquid component. The two interacting 
components are present at all temperatures in the underdoped regime, x < 0.22. 
At high temperatures T > T* the interaction leads to a modification 

of J, the effective exchange interaction in the spin liquid, while the Fermi liq- 
uid component remains intact. As the temperature is lowered, the growth of 
correlation length in the spin liquid results in stronger renormalization of the 
Fermi liquid component, and the formation of the Fermi arcs at temperatures 
T < T* ~ r™/3. The antiferromagnetic correlation length also becomes finite 
at low temperatures, which results in a gap for the spin excitation spectrum at 
Q = (tTjTt). The quantum critical point at x ~ 0.05 corresponds to a quantum 
transition between two different states of the spin liquid. 

We now return to the question of how close the measured spin liquid maps 
onto the Heisenberg spin liquid. While there is a one to one correspondence 
for the static spin susceptibility above T ~ T*, the difference becomes large at 
lower temperatures, where the correlation length ^ and the spin liquid gap l^sL 
saturate at a finite value. Presumably it is this difference that is responsible for 
the mild departure in concentration dependences deduced for f{x) and r™(x). 

4 Discussion 

Given the presence of both a spin liquid and a Fermi liquid in the normal state 
of underdoped superconductors, it is natural to consider their interaction as 
being responsible for their anomalous normal state properties and the transition 
to a gapped quantum state at T* and to superconductivity at T^- We begin 
by considering the extent to which the observed scaling behavior reflects that 
interaction, and what might be some of its other consequences. 

From an examination of the experimental results presented in the previous 
sections we draw the following conclusions: 

(i) The scaling with T™ of various transport properties of the non-Landau 
Fermi liquid tells us that the transport properties in the normal state of the 
cuprate superconductors are dominated by the magnetic coupling between the 
quasiparticles and the spin liquid excitations 

(ii) In the 1-2-3 materials, the scaling of T* and Ag, the d-wave energy 
gap found below T* for quasiparticles located in the "hot" regions of the Fermi 
surface, and below Tc for the remaining "cold" or nodal quasiparticles, tells us 
that the physical mechanism for their gaps and subsequent superconductivity 
must be magnetic. 

(iii) As may be seen in Fig. 43, the hot regions seen in the ARPES ex- 
periments are defined by their proximity to points on the Fermi surface that 
are separated by commensurate wavevectors, Q = (tTjTt); their physical loca- 
tion tells us that the magnetic interactions responsible for the gap must be 
peaked at Q. This is to be expected if the "glue" for the quasiparticle inter- 
action is provided by the spin liquid excitations of the Heisenberg model that 
exhibit the MMP form with a strong peak at Q, since their coupling to the 
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Fermi liquid yields an induced quasiparticle interaction that is likewise peaked 
at Q. Moreover, in agreement with the early models of spin-fluctuation-induced 
superconductivity[113, 114, 115] the resulting energy gaps will be d-wave. 

(iv) The pcnctration-dcpth measurements of the doping-dcpendcnce of the 
of the superconducting condensate tell us that where present in the underdoped 
materials, the spin liquid does not participate in their superconductivity. 

(v) The fact that the scaling relations below T"* for a pure Heiscnbcrg model 
are much better obeyed in the 2-1-4 materials than in the 1-2-3 materials tells 
us that the influence of the hot quasiparticles on the spin liquid is considerably 
weaker in the 2-1-4 materials. 

(vi) Since the spin liquid excitations in the two classes of materials are quite 
similar, these differences plausibly arise from band-structure effects[116, 117] 
that reduce the importance of hot spots on their Fermi surface. 

(vii) These last conclusions concerning the relative strengths of the coupling 
of hot quasiparticles to the spin liquid and their spin-fluctuation-induced quasi- 
particle interaction in the 1-2-3 and 2-1-4 materials provide a simple physical 
explanation of some of the major differences between their physical behavior of 
the two classes of materials. It is why the 2-1-4 materials show no indications of 
either a spin fluctuation gap or a quasiparticle gap opening up at T*, and may 
be related to their significantly lower superconducting transition temperatures. 

(viii) There are many indications that "stripy" behavior is more prevalent 
in the 2-1-4 materials; this may be due to the fact that the fluctuating domains 
become well established only when fairly strong AF correlations are present. For 
the 1-2-3 materials gaps open up at T ~ T™/3, so the maximum AF correlation 
length one expects to find is 3 lattice spacings; for the 2-1-4 materials, on the 
other hand, one finds AF correlation lengths ~ 8 lattice spacings near optimal 
doping and still longer correlation lengths are to be expected for the underdoped 
materials. 

(ix) Experiment indicates that at low temperatures T ~ T* ~ r'"/3, charge 
fluctuations, domains, and tendencies toward phase separation all play an im- 
portant role, although the two-component nature is present in the underdoped 
regime at all temperatures, and is therefore intrinsic. 

(x) Our two-fluid model provides a natural way to develop a phcnomenolog- 
ical description of stripe formation at low frequencies, since all that is required 
is to let f{x) vary periodically in space and time. 

We conclude this section with a brief discussion of some necessary elements 
of a microscopic description of the spin and Fermi liquids. We note first that for 
the spin liquid, the AF critical point at x ~ 0.05 is responsible for the observed 
z = 1 QC behavior. The spin liquid component at temperatures T > T* is well 
described to first approximation by the standard bosonic 0* Hertz-Millis theory 
of quantum critical bchavior[118, 119]. For ^ > 1 one is outside the critical 
region, in the vicinity of weak-coupling Gaussian fixed point, where fluctuation 
corrections are not important. The standard damping due to its interaction 
with fermions gives rise to the z = 2 mean field Ornstein-Zcrnike/MMP form. 
In the critical region 0* theory becomes equivalent to the QNLcr model[47, 48] 
(the two theories differ by an irrelevant operator). However, at low correlation 
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lengths the irrelevant operator of the cj)^ theory plays a crucial role, inducing 
a crossover to the mean field Gaussian regime[9] with ^ oc T~^/^ for N = oo 
version[120]. Such z = 1 to z = 2 mean field crossover cannot be obtained in 
the quantum non-linear cr-model with fermionic damping[121]. 

There appear to be two kinds of departures from the simple Heisenberg 
description. The first concerns the static bulk susceptibility that falls below the 
currently calculated values for the Heisenberg model for T < T™/3. Because this 
fall-off is essentially the same in both the 2-1-4 and 1-2-3 materials, it is likely not 
brought about by spin-fermion coupling, since that coupling is quite different 
for the two materials. It could represent an inadequacy in the calculations 
of the expected behavior of the Heisenberg model at these low temperatures. 
The second concerns the opening up of a gap in the low frequency spin liquid 
excitation spectrum at (tt, tt) that is seen as a minimum in TiT at 150K 
O.aiT™ in the NMR experiments on the 1-2-4 material and at 140^ ~ 0.23T" 
in the YBa2Cu3 06.63 material. It is plausible to conclude that it reflects a 
reduction in quasiparticlc damping brought about by the onset of the hot quasi- 
particle gapped behavior at T*. However, it would then seem that the transition 
to this new quantum state is not universal, i.e. does not always occur at a fixed 
fraction of T™. Indeed these NMR-determined onset temperatures, T*, scale, 
to first approximation with T^; thus for the 1-2-4 material one has T* ~ 2.1Tc, 
while for YBa2Cu306.63 one finds T* ~ 2.3Tc. Based on these estimates, and 
those obtained from the ARPES and STM experiments on 2212 members of the 
1-2-3 family , our candidate phase diagram for the 1-2-3 materials is given in 
Fig. 42b. It differs from the corresponding diagram for the 214 materials in 
the presence of a "high temperature" transition at T* to the gapped quantum 
state. 

In considering microscopic models for the evolution of Fermi liquid behavior 
with doping, we begin by recalling that high-Tc superconductors are doped 
charge-transfer insulators[122, 123]. Thus, an adequate starting point for a 
description of the Cu02 planes follows from their chemistry, and is the three- 
band Hubbard model composed of 3 dx2-y2 orbitals on the Cu site and 2 p^ y 
orbitals on the O sites[124, 125, 126]. The two key parameters of the three-band 
Hubbard model are are the charge- transfer energy Ect = £p — ^d, and the copper 
on-site Coulomb repulsion, Ud- In the limit of strong p-d hybridization and 
very strong on-site Cu Coulomb interaction, the model gives rise to 5 different 
bands: an antibonding band, further split into upper Hubbard band (Cu d^" 
configuration), and lower Hubbard band (Cu dP configuration) by the Coulomb 
interaction, a bonding band, further split into Zhang- Rice (ZR)[35] singlet and 
triplet bands, and the non-bonding oxygen band. 

Following Anderson[97], most theoretical studies have considered a simpler 
one band Hubbard model, which has fewer parameters and assumes that the 
low-energy physics of the electronic states in the vicinity of the charge-transfer 
gap is described well by taking a relevant low-energy subset of this complex 
band structure, the ZR singlet band and the upper Hubbard band, separated 
by [/ ~ Ect ^ l-7eV of the order of the charge-transfer gap. Numerical studies 
of the three-band Hubbard model[127] support such reduction. The two low- 
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energy bands of the effective single-band Hubbard model are strongly correlated 
hybridized copper-oxygen bands, with the upper Hubbard band having mostly 
oxygen character[128, 129, 116, 117]. 

However, on the basis of numerical LDA calculations, it has been suggested 
that a single-band tight binding model may not provide a sufficiently good quali- 
tative description of the many bands involved[116, 117]. Moreover, we have seen 
that Xdd contains the spin liquid component localized on the copper sites, while 
Xpd and Xpp pick up contributions from the holes on the oxygen sites. The differ- 
ence in relevant form factors could potentially explain many anomalies observed 
for NMR relaxation rates [36] and inelastic neutron scattering experiments, as 
we have discussed in the previous section. 

The analogy with two-fluid model for heavy fermions[101, 102] can be better 
understood from an alternative approach to three-band Hubbard model[124, 
125, 126], that is perturbative in p-d hybridization, tpd- In case of exactly one 
hole per Cu and strong Hubbard Coulomb on-site repulsion Ud, the low-energy 
Hamiltonian for the three-band Hubbard model reduces to the Heisenberg spin 
Hamiltonian with exchange coupling 



since the charge excitations have a gap E^j, ~ 1.7eV. Numerical calculations for 
the three-band Hubbard model[127] extract the value Jh — O-lSeV, consistent 
with that observed in neutron scattering experiments. The oxygen holes and 
the nearest Cu spin interact via a strong Kondo exchange interaction, which 
can be obtained in perturbation theory in tpd/Ecr, 



The perturbative analysis of the three-band Hubbard model reduces it to the 
spin-fcrmion model via the Schricffer-Wolff[130] transformation. 



lation operators, respectively, for holes in oxygen p bands. The Kondo inter- 
action is, in general, non-local[125]. The Hamiltonian Eq.(112) is the Kondo- 
Heisenberg model, known for the periodic Anderson model in many physical 
systems, including magnetic semiconductors[131], heavy fermions (the Kondo 
lattice model) [130], and high-Tc superconductors[124, 125]. The heavy fermion 
materials and high-Tc materials correspond to different limits of the model 
Eq.(112). In heavy fermion superconductors the Kondo interaction Jx is ini- 
tially comparatively small, and gets renormalized by the conduction band as the 
temperature is lowered, leading to the onset of coherence at T* , and the two- fluid 
heavy fermion behavior[101, 102, 132] with a temperature-dependent fraction 




(110) 




(111) 
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of spin liquid (spin impurity contribution). The high-Tc superconductors, on 
the other hand, are in the strong-coupHng Zhang-Rice [35] Hmit of the same mi- 
croscopic model, Eq.(112). The Kondo interaction is very large, Jk ^ O.TeV, 
which leads to a formation and propagation of the "Zhang-Rice singlets", re- 
ducing the microscopic Hamiltonian to an effective one-band t — J model [97]. 
The projection operators of the Zhang-Rice singlet [35] automatically imply that 
Xpd = Xpp = 0- Thus, as suggested by MMP[40], different nuclei then feel the 
presence of a single temperature-dependent spin component residing on copper, 
or Xdd- The Zhang-Rice singlet picture (and the t-J model), however, inevitably 
gets violated at higher doping, as a single large Fermi surface develops, that cor- 
responds to antibonding p-d hybridized band. A linear decrease of effective J{x) 
with increased doping x is "high-energy physics" , as is evident from, for exam- 
ple, an adiabatic treatment of the t-J model[133]. Thus, the fraction of spin 
liquid [3] observed in "low energy" experiments on the cuprate superconductors 
will be temperature-independent. In heavy fermion metals it is the mixed c — / 
contribution, Xcf{T) = (1 — f{T))xFL{T), that gives rise to the heavy Fermi 
liquid contribution[102]. In high-temperature superconductors, on the other 
hand, one normally uses the Zhang-Rice singlet band and the upper Hubbard 
band of the three-band Hubbard model as lower and upper Hubbard bands of 
the effective single-band Hubbard Hamiltonian. Using the analogy the heavy 
fermion two-fluid model[101, 102] and the two-fluid model for high-T^ super- 
conductors, the Fermi liquid contribution in the cuprates should arise from Xpd 
and Xpp parts of the bulk spin susceptibility, which, strictly speaking, indicates 
a violation of the Zhang-Rice singlet picture at high doping, given by Eqs. (50), 
(51). Since the effective mass of carriers is not heavy in high-Tc superconduc- 
tors, one would expect both Xpp and Xpd contributions to be present. This 
result contrasts strongly with the pure Zhang-Rice [35] limit in which all contri- 
butions arise from Xrfd[40]. The key answer to this very important question is 
held by the new Knight shift experiments[102, 39], since then Knight shifts on 
different nuclei should see a different temperature dependence, at least below 
the superconducting temperature, Tc- 



5 Conclusion 

Our review of scaling behavior in the cuprate superconductors makes evident the 
critical role that the localization of Cu spins by strong on-site repulsive Coulomb 
interactions plays in determining the properties of the magnetically underdoped 
superconductors, those whose planar hole concentrations are < 0.20. While 
there is much more in the way of details that experiment can and will provide, 
from additional c-axis NMR experiments on 1-2-3 materials that enable one to 
follow the emergence of a distinct Fermi liquid contribution, to ARPES, NMR 
, and STM experiments that fill in the details of the emergence of the gapped 
hot quasiparticle state at T* , we believe that the overall low-energy physical 
picture of a normal state described by two fluids whose interaction leads to its 
remarkable properties has now been firmly established. 
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It is equally clear that much much more remains to be done, beginning with 
the development of a better microscopic foundation for the physical picture 
that experiment has "forced" upon us, and an improved treatment of transport 
processes based upon it. We mention here but a few of the many remaining 
challenges and mysteries that surround these fascinating materials: 

(i) There are striking similarities in the doping dependence of Tc in the 1-2-3 
and 2-1-4 materials despite the difference in their spin-fermion couplings and 
the magnitudes of their respective transition temperatures. Does this tell us 
that these differences are primarily in how the hot quasiparticles couple to spin 
fluctuations and that the spin-fermion couplings for the cold quasiparticles that 
bring on their superconductivity can somehow be scaled for these two classes of 
materials? 

(ii) Although the spin liquid ceases to be a significant "player" for doping 
levels greater than, say, 0.20, the existence of a smooth cross-over in most phys- 
ical properties of the normal state, including the superconducting transition, as 
one goes from magnetically underdopcd to magnetically ovcrdoped materials, 
suggests that its disappearance must be accompanied by the emergence of suffi- 
ciently strong antifcrromagnetic correlations in the Fermi liquid that enable its 
properties to merge smoothly with those induced by the spin liquid. Can these 
be detected experimentally and understood microscopically? Might the disap- 
pearance of the spin liquid be accompanied by the emergence of a single-band 
Hubbard model that is valid for the overdoped materials? 

(iii) Do the gapped hot quasiparticles found below T* in the 1-2-3 materials 
become effectively localized and so become part of the spin liquid? Is there any 
way to make their gapping more productive, in terms of their becoming part of 
the superconducting condensate, thereby opening the way to a superconducting 
transition at the considerably higher transition temperature, T*? 

(iv) Can the lessons learned from the cupratcs be transferred to other ma- 
terials, and so lead over time to the discovery of still higher temperature super- 
conductors? 
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